General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



Cft. 933fcl 


0,0*11 A0l>e TO Ttf£ ft>ISUC 


STD RESEARCH CORPORATION 

BOX 4-1 27, CATALINA STATION * PASADENA, CALIFORNIA 91106 

TELEPHONE: (213) 684-1771 



Contract No. NAS 2-3580 


STD- 68-2 


NOVEL METHOD FOR DETERMINATION OF ENERGY -LOSS FACTORS FOR 
SLOW ELECTRONS IN HOT GASES 


INTERIM REPORT 

Covering the period 28 April 1966-30 June 1968 


N 

s 

z 

K 

O 


N69-3541 3 



(THRU) 

J 

(Cf DEI 
(CATEOORY) 


% 


Prepared for 

National Aeronautics and Space Administration 
Ames Research Center 
Magnetoplasmadynamics Branch 
(Howard Stine, Technical Monitor, Mail Station N-229-2) 
Moffett Field, California 94035 


1 August 1968 


f RECEDING PAGE clank KOT RL/wLO. 


CONTENTS 


SUMMARY 

INTRODUCTION 

ELECTRON ENERGY TRANSPORT 

AND THE ENERGY -LOSS FACTOR 

DESIGN OF THE EXPERIMENT 

Physical Constraints Placed on the Experiment 

Operating Conditions 

Experimental Apparatus 

ENERGY STATE OF THE GAS 

EXPERIMENTAL RESULTS 

Experiments with Argon 

Experiments with Nitrogen 

Experiments with Nitrogen- Oxygen Mixtures . 

DISCUSSION OF RESULTS AND CONCLUSIONS . . 

ACKNOWLEDGMENTS 

APPENDICES 

A. Dimensions of Waveguide and 

Estimate of Wave Impedance 

B. Sensitivity of the Measurements .... 

C. Experimental Probe Characteristics; 

Surveys of Electron r.nd Ion Currents 
and Electron Temperatures; Typical 
Rotational and Vibrational Spectra . . . 

REFERENCES 

TABLES 


Page 

1 

2 

3 

11 

15 

25 

29 

38 

45 

46 
49 

52 

53 

57 
53 

58 
64 


67 

71 

76 


FIGURES . 


85 



CEDING PAGE 


Table Title Page 

I. Summary of Constraints . . . . . 7 6 


II. Typical Plasma Temperatures in °K 

for Arc Heated Pure Nitrogen Expanded 
to a Pressure of Approximately 400 


Microns 78 

III. Typical Energy-Loss Factors for 

Argon at a Gas Temperature 

Between 100 and 200 °K 79 

IV. Typical Electron Temperature 

Measurements as a Function of 

rf Power for Nitrogen 80 

V. Typical Energy-Loss Factors for 

Nitrogen as a Function of Gas 

Temperature 81 


VI. Typical Energy-Loss Factors for 

Oxygen-Nitrogen Mixtures as a 
Function of Composition and Gas 

and Electron Temperatures 82 

VII. Dependence of Energy Loss 

Factor for Pure Oxygen, Derived 
from Mixture Data by Mixture 
Rule* and Assuming that only N£ 
and O 2 Are Present in the Mixture, 

on the Oxygen Content of the Mixture 84 


v 



fiilCEDiKG PAGE CLAN- KOT PiUv^p. 



LIST OF FIGURES 


Figure 

Title 

Page 

i. 

Schematic of microwave plasma 
illumination system 

85 

2. 

Schematic of special waveguide for 

plasma illumination 

86 

3. 

Energy-loss factor for N 2 -O 2 mixtures 

at T * T .. * T _ A « 1700 °K as a 
g vib e , E = 0 

function of oxygen mole fraction 

87 

A- 1 . 

Cross section of the waveguide 

88 

C- 1. 

Langmuir probe characteristics for 
nitrogen plasmajet with 2. 4 watts and 
0 rf power (traced twice) 

89 

C-2. 

Langmuir probe characteristics for 
nitrogen plasmajet with 2, 4 watts and 
0 watts rf (traced five times) 

90 

C-3. 

Radial survey of electron and ion 
current to Langmuir probe obtained by 
continuous traverses of probe through 
plasmajet axis and probe characteristics 
on axis 

91 

C-4. 

Axial survey of electron and ion current 
to Langmuir probe obtained by continuous 
traverses of probe along plasmajet axis . . , 

92 

C-5. 

Radial distribution of T 

e 

93 

C-6. 

Langmuir probe characteristics along 
axis of plasmajet obtained at t u 
intervals from nozzle exit to 20" 
downstream 

94 


vii 




LIST OF FIGURES (Concluded) 


Figure Title Page 

C-7. Electron temperature distribution 

in waveguide (axial) 95 

C-8. Effect of rf on probe electron current 

on jet axis (axial traverse) 96 

C-9. Rotational band of N^" (0, 0) used to 

determine rotational temperature 97 

C-10, Band profiles of (0, i) progression 

at T vib * 2000 ° K > T vib * T rot used 

to determine vibrational temperature 98 


viii 



LIST OF SYMBOLS 


a Width of waveguide 

Cross-sectional area of waveguide 

b Height of waveguide 

B Magnetic induction 


d Thickness of plasma, equation (A-3) 

Dp Diameter of Langmuir probe 


e Charge of the electron (a positive number as defined), 

e = - e 
e 

E Applied (rf) electric field, E = electric field vector 

E q Amplitude of applied (rf) electric field 

AE q Experimental error in E q , Appendix B 

Ep Characteristic plasma field, equation (31) 


f 

f 

(3f e /at) c 

f 

p 


®e k 
^eK^ o 


Frequency (cycles per second) 

Distribution function for electrons 

Rate of change of f due to collisions 

e 1/2 

Critical plasma frequency (Hz), f = 9 n & ' 

Relative speed between an electron and a molecule of 
species k 

Relative speed g averaged over a Maxwellian 
distribution 


J 

e 

k 

k 

o 


Total electron current 

Boltzmann 1 s constant or (Appendix A) wavenumber, 2tt/\ 
Wavenumber, equation (A-l) of Appendix A 


m 

e 

m 

K 


Mass of the electron 

Mass of a molecule of species k 

Mole fraction of oxygen in experiments with N^-O^ 
mixtures 


IX 



n 

n 


a 


e 


n 


K 


Neutral particle density 
Number density of electrons 
Number density of species k 


p e 

p 

(AP) 

P 


min 


nm 


Pressure of the electron gas 

Applied rf power (equation (19) and Appendix A) 

Minimum power required, defined following equation (29) 

Power transmitted in propagating mode of waveguide, 
equation (A- 5) of Appendix A 



Q 


e* 


Q 


e k, s 


Energy flux vector for electrons, equation (1) 

Collision cross-section for momentum transfer between 
electrons and neutral particles of species a 

Effective cross-section for momentum transfer between 
electrons and species k , equation (14) 

Effective cross-section for energy transfer between 
electrons and the s degree of freedom (translational, 
rotational or vibrational) of species k , equation (12) 


P< 2 > 


[r (2) ] , 

e J el 

[R (2) ]. , 

e mel 


R 

R 


M 


Collisional energy exchange integral, equation (2) and 
equation (4) 

Collisional energy exchange integral, elastic collisions 
only 

Collisional energy exchange integral, inelastic collisions 
only 

Ratio defined by equation (21) 

Radiative energy loss rate (ergs/sec - cm ) 


s. 

i 


Rate of reaction (number/ sec - 


cm 3 ), 


equation (4) 


t 

T 

T e 

AT e 

< AT e>min 


T 

e, E = 0 
T 

K 


Time 

Temperature 

Electron temperature; T^ denotes an initial value 
Increment of electron temperature, AT = T - T , 

v C “A 

Minimum permissible change of electron temperature, 
equation (29) 

Electron temperature with no applied electric field 
Temperature of species k ; T = effective gas temperature 


x 



k trans 
L rot 

n 

k vib 


elec 


k, s 


T n 

T? 

J 


Temperature associated with translational motion, T 
Temperature of rotational degrees of freedom 
Temperature of vibrational degrees of freedom 
Stagnation temperature 

Temperature of electronic degrees of freedom 

T k, rot ' T K,vib- or T «,tr < 8ee equation (8) ) 
Reference temperature, equation (9) 

Excitation temperature for a class j of energy states 
(electronic, vibrational, rotational) 


tr 


u 


U 


Particle velocity for electrons 
Mass average velocity of gas 


w 


W 


Particle velocity (drift) relative to mass average velocity 
Average of w over velocity distribution, W = (w ) 

G G C 


X 

x 

AX 


Coordinate along direction of flow 


X = T 


Appendix B 


Experimental error in X , Appendix B 


y. z 


Z 

Z, 

k, n, m 


Coordinates normal to direction of flow 

Impedance, equations (19) and (A-7), Z = (|i/e)*'^ 
Defined in equation (A- 2) 


6 

A 5 

6 

K 

6 1 
k, el 



Electron energy loss factor 
Error in 6 , Appendix B 

Energy loss factor, average fraction of energy lost per 
collision between an electron and a molecule of species k 

Energy loss factor for elastic collisions, 6 , = 2 m /m 

K, GI G K 

Energy loss factor for the s class of energy states 
(translational, rotational or vibrational) for species k ; 
see equation (13) 


Energy loss factor, electron-electron collisions 
Effective energy loss factor for a gas mixture 

5 H 6 k,s • < T e;X" T e- T K,. )/(T e- T n> 


K 8 


XI 



e 

€ . 
1 

1 \ 


0 

0 

0 


V 

rot 

elec 


Quantity appearing in equation (17) 
Energy per reaction, equation (4) 

r| - (1 - */ cj 2 ) 1 7 ^ , Appendix A 


Charactei istic temperature for vibrational excitation 
Characteristic temperature for rotational excitation 
Characteristic temperature for electronic excitation 


\ Wavelength, Appendix A 

\ Debye length, equation (26) 


v eff 


V 


te 


v 

v 


t 

e, 


K 


V 

e, k, s 

n 


Effective number of collisions per second for electrons 

Total effective collision frequency for energy transfer 
(electrons) 

Defined in equation (6) 

Collision frequency for electrons with molecules of 
species k 

Defined by equation (12) 

Quantity appearing in equation (17) 

Stress tensor for electrons, t t = p (w w ) 

e e ' e e ' 


P 


Reflectivity, equation (A-3) 


T 

e, e 

T 

e, e, e 

T 

e, k 

T 

e, k, s 


u> 



Effective collision time for momentum transfer 
among electrons 

Effective collision time for energy transfer among 
electrons 

Effective collision time for momentum transfer between 

electrons and particles of species k, (4/3)t = l/v 

e, k e, i 

Defined by equation (12) 


* = ft + Za.e. 

i 1 1 

Angular frequency of applied electric field 
Plasma frequency, Up = 2irfp = 2 tt fb 
Ohmic heating, equation (38) 


) 


xii 



Subscript 


( ) Quantities associated with electrons 

' e 

( ) , Values for elastic collisions 

' 'ei 

( ). . Values for inelastic collisions 

'mel 

( ) Average over velocity distribution 

( ) Average for a Maxwellian distribution 

( ) Quantities associated with the s class of energy states; 

8 s identifies translational, rotational or vibrational 

degrees of freedom 

( ) Quantities associated with species k 


xm 



NOVEL METHOD FOR DETERMINATION OF ENERGY -LOSS FACTORS FOR 

SLOW ELECTRONS IN HOT GASES 

by S. T, Demetria.de s 

STD Research Corporation, Pasadena, California 91106 


SUMMARY 


The energy-loss factor for slow electrons in hot gases is determined 
by a new method that employs a high-frequency (2. 45 GHz) electric field to 
elevate the temperature of the electrons above the temperature of the gas 
and a Langmuir probe to determine the electron temperature. The electron 
energy-loss factor 8^ in a given gas k is then determined by measuring 
the rate of change of the electron temperature with the high-frequency power 
used to illuminate the plasma. The determination of 8^ by this method 
does not require knowledge of the collision cross-sections for momentum 
transfer or of other transport coefficients or of the number densities of the 
electrons and neutral particles present. An arc heater was used to generate 
the hot gas. The energy state of the gas was determined spectroscopically. 


The values of 8^ obtained by this method are in agreement with 

theory for monatomic gases (e. g., argon). The values of the energy-loss 

factors 8^ obtained fcr nitrogen and 5^^ obtained for N^-O^ m * xtures 

(with oxygen concentrations up to 50% ) are a reasonable extrapolation of 

the room temperature data available in the literature. It was found that 

- 4 o 

8^ ranges approximately from 3X10 at a gas temperature of 1700 K 

to 7 X 10' 4 at 5000 °K. The determination of electron energy-loss factors 

in mixtures is further complicated by the large number of variables; 

nevertheless for electrons in the temperature range between 2000 °K and 

2500 °K in air at temperatures in the vicinity of 1700 °K it was found that 

5 . *3X 10" 3 . 

air 


It is concluded that this is a valid and accurate method for deter- 
mining energy-loss factors since the experimentally obtained values are in 
essential agreement with theoretical predictions for various gases over a 
range of more than two orders of magnitude in the loss factor. 



INTRODUCTION 


Many non-linear phenomena in plasma dynamics and electro- 
magnetism as well as the processes of excitation or de -excitation of 
atomic or molecular particles by electron impact are characterized, and 
to a large extent, controlled, by the electron energy-loss factor, a quantity 
that determines the rate of energy transfer from electrons to heavy parti- 
cles and expresses the fraction of energy that an electron loses, on the 
average, per collision with the heavy particles. 

Development of a simple, reliable and direct method for the 
determination of electron energy-loss factors in hot gases that permits 
presentation of the results in a concise and convenient form (e. g. , ai a 
function of electron and gas temperature) is the purpose of this investi- 
gation. 

The energy-loss factor for slow electrons in hot gases is determined 
by a new method that employs a high-frequency electric field to elevate 
the temperature of the electrons above the temperature of the gas and a 
Langmuir probe to determine the electron temperature. The electron 
energy-loss factor 6^ in a given gas k is then determined by measuring 
the rate of change of the electron temperature with the high-frequency 
power used to illuminate the plasma (ref. 1). 

The experimental results obtained so far offer evidence that this 
method is capable of meeting the objectives: it is demonstrated that the 
proposed method works and gives results that can withstand all reasonable 
checks for validity, that the energy state of the gas can be adequately 
established and that the results can be obtained and presented in a concise 
and convenient form. 

In this method for determining the electron energy-loss factor in 
hot gases all important measurements are dc and the energy-loss factor 
is uncoupled from (a) the collision frequency for momentum transfer, or 
(b) the drift velocity, the mobility and the electron diffusion coefficient 
and hence can be determined independently of particle densities, collision 
cross-sections for momentum transfer, and ether transport parameters. 
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Experimental data on the value of 6^ as a function of electron and 
gas temperature obtained by this technique are in a form consistent with 
the powerful new theoretical treatments of non-linear plasma dynamics now 
becoming available (e. g. , such as the theoretical work of Argyropoulos 
et al. , ref. 2) and are essential for obtaining accurate quantitative results 
from these theories. 

For these reasons a recapitulation of some of the pertinent basic 
theoretical expressions of non-linear plasma dynamics is appropriate and 
useful. 


ELECTRON ENERGY TRANSPORT 
AND THE ENERGY -LOSS FACTOR 


Success in any scheme to describe quantitatively the electron- 
energy dependent non-linearities in a plasma rests on the degree of under- 
standing that has been achieved concerning the mechanisms of electron 
energy transport in high-temperature gases. 

A consistent scheme for obtaining a closed set of transport equations 

in a plasma was given by Demetriades and Argyropoulos (ref. 3). Using 

their notation, repeated here in the list of symbols, it can be shown 

(ref. 4) that the third moment of the Boltzmann equation for electrons 

(i. e. , the equation obtained by multiplying the Boltzmann equation by 
i 2 

and integrating over velocity space) results in the so-called 
"thermal energy balance" equation which is given by 


a 

St 


<2 P. 


p w 

r e e 


) + V • (-| p e U + ~q e ) + tr e : 

. (E + u X B ) - 

m ' ' dt 

l e 


VU = 

+ r< 2 > 

e 


(i) 


Equation (1) is also obtained by subtracting the dot product of the 
momentum equation for the electrons and the flow velocity U from the 
total energy equation for the electrons (whence the term "thermal energy 
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balance equation). In equation (l)the subscript e refers to electrons 

( 2 ) 

and the collisional energy exchange integral ' is defined by the general 
expression* 

„(2) _ C i 2 , 9f e, 

R e = ) 7 m e w e dw e < 2 > 

w 

e 

where f is the electron distribution function, 
e 

The elastic part of this integral (i. e. , the translational energy 
exchange) can be readily obtained by the methods of reference 3 or refer- 
ences 4, 5 and 6 and reduces exactly to the form 

[R (2) ] , = - l kn y (2m /m )(T - T )t -1 (3) 

1 e J el 2 e v e / k' v e k' e, k v 7 

under the assumption that the isotropic part of the distribution functions 
both for the electrons and for all heavy components k is Maxwellian. Thus 
equation (3) as written above is also a result of the general elastic collision 
integral expression for ; given by equation (3) of reference 3. 

In addition to the elastic part given by equation (3), one generally 
considers the following inelastic collisional processes of energy exchange 
between free electrons and heavy particles: 

(1) Excitation of rotational and vibrational levels in 
molecular species 

(2) Excitation of electronic states of the heavy particles 

(3) Collisions leading to dissociation, ionization cr recombination 
or the creation cf metastable species with lifetimes much 
longer than the lifetimes of the excited particles of item (2) 

The identity of this definition and that appearing as equation (3) of 
reference 3 may be established by use of known properties of molecular 
collisions. See, for example, Hirschfelder, J, O, ; Curtiss, C. F. ; and 
Bird, R. B. : Molecular Theory of Gases and Liquids, J. Wiley Sons, 

New York, (1954) p. 460. 
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(4) "Second order" impacts that result in the transfer of the 

energy of the excited or metastable particle to the incoming 
electron 

As shown by the authors of reference 7, the inelastic energy losses 
to rotational and vibrational levels can be adequately taken into account by 
substituting 6^ for 2m e /m^ in equation (3) and using a measured value 
of 6^ (such measurements are described in reference 1), 

Concerning the excitation of electronic states, use can be made of 
the fact that in most plasmas of interest to this investigation the excited 
levels of atoms or molecules are generally in instantaneous equilibrium 
with the free electrons, [in these plasmas, while the electrons, atoms and 
molecules experience energy exchange among one another, the ultimate loss 
of energy from the system is due to radiation. ] Therefore, the appropriate 
portion of the collisional energy losses of free electrons (i. e. , that portion 
of the energy lost by an electron upon impact with an atom or molecule that 
goes into exciting electronic levels of the heavy particle) can be equated 
to the radiative losses R from excited states of the heavy particles at the 
electron temperature (see for example reference 8). 

Finally, collisional energy losses due to ionizing or dissociating 
collisions or collisions leading to the creation of metastable states (free 
radicals, etc. ) or energy gains due to the reverse reaction mechanisms, 
are expressed in terms of the rate s^ of each reaction (in units of number 
of reactions per unit volume per unit time) and each characteristic energy 
or potential (in units of energy per reaction). The sign of depends 
of course on whether the electron loses or gains energy by that particular 
mechanism 

Thus the total (elastic + inelastic) collisional exchange of energy 
between electrons and heavy particles can be written as the sum of equation 
(3) and the inelastic contributions just enumerated, in the form 


R 


( 2 ) 





2 

m w 


e e 


w 


e 


0f _ 

<-5 f> c d *e 


5 



kn Yj & (T - T ) t -1 
e u k e k e, k 

t/ ' 


l s.e. 

V 11 


The quantity 6^ accounts for the energy lost by the electron in elastic 

collisions (i. e, , involving translational degrees of freedom only) as well 

as in some inelastic collisions (i. e. , involving rotational and vibrational 

degrees of freedom of the heavy particles). The quantity 6^ is thus a 

characteristic energy-loss parameter that has the meaning of an average 

3 

fraction of the energy difference 7 k (T - T ) transferred in collisions 

between the electron and the heavy particles k within a time T e K • For 

example, when the heavy particles are at equilibrium and have only one , 

well-defined temperature T^ = T^ , the difference between the average 

internal energy of the electrons and the average internal energy of the 

3 

heavy component k is ^ k (T^ - T ) . Then if the heavy particles are, 

for example, monatomic, optically unexcited and unexcitable, chemically 

stable and inert (i. e. , if they behave like hard inert spheres) the electron 

impact is perfectly elastic and the elastic energy loss factor ^ K e i the 

fraction of this energy difference that is lost on the average by an electron 

in each elastic collision with component k . Thus 6 , = Zm /m . 

" ■ 1,1 1 ■ k, cl e k 

In a single -component gas or plasma at equilibrium where again 

and Tg are equal, and R as well as Ss.c. are known, the appropriate 

inelastic energy-loss factor 6^ can be determined experimentally, for 

example, by obtaining all other quantities and using equations {1) and (4). 


In a multicomponent gas or plasma where all the heavy components 

are at equilibrium at the same temperature T = T , we can write 

k § 

equation (4) in the form 


kn v.6 ,, (T - T ) - R - Y. s.e. 
e t eff ' e g 1 v l l 


where v. is defined in reference 3 as 


»t ■ ? 


• 4 l 


6 



and 6 


is defined in reference 1 as 


eff 


6 


eff 


<l/v t > l 6 k 

K 


-1 

T 

e, k 


(7) 


Again 6^ can be obtained experimentally [either directly from equations 
(1) and (5) by obtaining all other quantities or from the 6 *s for each 

K 

component and the mixture rule given by equation (7)]. 


Finally, in a multicomponent gas or plasma where the heavy com- 
ponents are at different translational, rotational and vibrational tempera- 
tures T , T and T respectively, we can write 

k ^ tr k j rot K, VXD 


R 


( 2 ) 



E E 6 k, s < T e 

K S 5 


- T ) T 
K , s e, 


« 

R 



(8) 


where the subscript s identifies translational (tr>, rotational (rot) or 
vibrational (vib) temperatures or corresponding contributions to 8 . 

Once more the values of the coefficients 6 can be obtained 
. K, S 

experimentally provided R and 2 s.e^ are known or measurable. An 
effective energy-loss factor can now be defined as 


6 


eff 


E E * 

K S 


-i 

T 

e, k 



(9) 


where the double summation requires two normalization factors (i. e. , 
and T g - T ) and T n is any convenient reference temperature. Then 
we obtain 


R 


(2) 


3 

1 


k n 


v t 6 eff 


< T e " 


T ) - 
n' 


R 


E 


S.€ . 
1 X 


(10) 


It should be remembered that the first term on the right-hand side of 
equation (10) represents the energy exchange between electrons and heavy 
particles in collisions that excite translational, vibrational and rotational 
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energy states. Hence this term involves all of the elastic and part of the 
inelastic energy exchange. 

Note that the summation over the classes s of energy states where 
s = tr, rot, vib is equivalent to treating each class s of energy states of 
each species as a distinct and separate component. When this is done, 
the first term on the right-hand side of equation (8) is analogous to and a 
direct extension of the simple expression for the translational energy 
transport, equation (3) above, derived from the theory presented by 
Demetriades and Argyropoulos, reference 3. The assignment of a tempera- 
ture T to each state of each plasma component and the treatment of 

K, S 

each state of each component as a different species is permitted on the basis 
of energy considerations. On the other hand, the collision frequency T e ” 
is excluded from summation over energy states since it characterizes 
momentum transfer only. 


As an example, let us consider a mixture of nitrogen and oxygen 

where Tkt_ + _ = T^t . = T n „ , = T 0 ^ . while T-, 

iN 2 , tr iN 2 , rot O 2 , tr C>2 , rot N 2 , vib O 2 , vib 

T ., and 6, r . * , < 6^ T , < 6^. , « 6» T < 6 ., . 

vib N 2 , tr O 2 , tr N 2 , rot O 2 , rot N 2 , vib O 2 , vib 

Then by equation (9) the energy-loss factor for the mixture becomes 


<• „ e vib 1 . V , -1 

°eff ~ T - T v, L-t °k, vib e, k 

e n t k ’ 

where k = N_ , 0 0 . Then if we choose T = T for normalization 
2 ' 2 n vib 

purposes, we obtain expressions identical to equations (5) and (7) where, 
in this case, the appropriate heavy particle temperature, T^ , is the 
vibrational temperature. 


It is important to note at this point that 
known species) is much smaller than unity. 


(and 6 for all 

' K 


,c By class of energy states we mean the energy states that are associated 
with a particular degree of freedom. For example the vibrational energy 
class contains all the energy states that are associated with the vibrational 
degrees of freedom. 
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We can nuw proceed to make some statements illustrating and 
clarifying the use of equations (1) and (5), (8) or (10). 

Use of the term R in equations (5), (8) or (10) implies that the 
quantity 6 , or more generally 6 excludes collisional excitation of 

K K , S 

9 t 

bound electrons by free electrons. Use of the term R is based on the 
assumption that the population of the bound electronic levels is at equili- 
brium at the temperature of the free electrons T g or, in other words, 
that the energy lost by the free electrons in the process of excitation of 
optical levels of heavy particles is equal to the energy lost by these heavy 
particles through radiation. Clearly, the quantity R is a complicated 
function of the geometry of the plasma (refs. 9 and 10) and the relative 
importance of the term R compared to the term containing the energy loss 
factor is greater in monatomic gases (where 6^ is close to the elastic 
value 2m /m and therefore quite small) than in diatomic gases (where 
6 is usually many times greater than 6 1 ). Usually, however, even 

for monatomic gases there are ranges of T g and where the term 
containing 6^ in equations (5), (8) or (10) is dominant and R can be 
safely neglected (e. g. , see reference 11). 

The other important term, S s.e. , appearing in equations (5), 

ill 

(8) or (10) may be the source of considerable confusion if it is not properly 
taken into account. For example it is easily seen that it can be both a 
source and a sink of electron energy (as would be the case for example when 
free electrons are picking up energy by becoming collision partners with 
metastables or third bodies in recombination reactions while at the same 
time they become collision partners in dissociative reactions). A reason- 
able estimate of the relative magnitude of this term within the control 
volume of the experiment will always be necessary when use is made of 
these equations. 

A final point that must be clarified in this section is the relationship 
between collision cross-sections for momentum transfer and collision 
cross-sections for energy transfer. We define an "effective collision 
frequency for energy transfer" between electrons and the s-class oi energy 
states of heavy particles of species k in a manner completely analogous 
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to the definition in reference 3 of an "effective collision frequency for 
momentum transfer" t , by the expression 

G j P 


-1 

T 

e, k, s 



^eK^o 


Q 


e, k, s 


4 

T V 

e, k, s 


( 12 ) 


where (g eK ) Q is the mean approach or relative speed between the electron 
and component k (for Maxwellian velocity distribution of the two com- 
ponents at their respective kinetic or translational temperatures) and 
Q is the "effective collision cross-section for energy transfer" 

between electrons and the s-class of energy states of species k . This 
definition is entirely consistent with the transport theory formulation and 
methodology of reference 3. 

( 2 ) 

Now we consider the meaning of the term R^ ' . It can be thought 

of as a collisional energy source or sink for the electrons (in the same 

manner for example, that which is treated in reference 3, can be 

thought of as a collisional force* ) and represents a power density. 

Therefore, in a manner analogous to the make-up of R^ , the first term 

on the right-hand side of equation (10) must be made up of an energy 

difference term, a number density term, and a collision frequency. A 

( 2 ) 

look at the elastic expression for R' , equation (3), should suffice to 

6 3 

con/ince us that the energy difference term is k(T - T ) , the number 
density term is n g and therefore the effective collision frequency for 
energy transfer can also be defined by 


-1 

T 

e, k, s 


6 

K, 


S 


-1 

T 

e, k 


4 

I 


k, s e, K 


(13) 


Therefore, using the result of reference 3 for the effective collision 
frequency for momentum transfer given by t = rr n (g ) Q 

v | K j K C K 

we obtain 

* R^ is defined by an integral similar to that for R^ , equation (2), 

e ^ 12 e 
except that m w replaces 7 m w . See equation (3) of reference 3. 

© 0 L 6 6 
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( 14 ) 


e, k, a 


From equations (10) and (13) we see that the total effective collision 

frequency for energy transfer v tg is equal to v t and not to the simple 

sum S 2 t * unless all T = T . Since electron collisions with 
k s e, k, s k, s n 

each component k in the energy class s contribute a fraction of the total 
energy exchange that must be weighted by the temperature ratio 
(T - T )/ (T - T ) , it follows that the total effective collision 
frequency for energy transfer is given by the expression 


v 


te 


zz 


-1 

T 

e, k, s 



v t 6 eff 


(15) 


It must be emphasized that equations (14) and (15) are really 
consequences of the definitions given by equations (12) and (13) and although 
they add nothing new to our knowledge they help maintain the consistency of 
our results and underscore the points of similarity and/or difference with 
other theories or investigations 


DESIGN OF THE EXPERIMENT 

Under steady-state conditions (3p e /3t = 0 , 3U/3t = 0 ) in a flow 

field with negligible gradients of axial velocity, electron concentration and 
temperature (i. e. , V • [3p U/2 + q j +ir : VU = 3 U • V kn T /2 + 

6 6 6 v 

3pV*IJ/2 + V‘"q' + p V - IJ - : VU « 0 since 3n /3x = 3T / 3x 

1 e ' ^e r e e e e 

= 3U x /3x ~ 0) the electron thermal energy balance, equation (1), becomes 

T 0 • (IT + U X B) = | * n e v t 6 eff (T e - T n ) + R + Z (“I 

where 6^ is defined b’ r equation (9). When U X B is negligible and the 
electric field is constant with time we obtain, in the first approximation 
(ref. 3) 


11 



(16a) 


E = 


2 r 2 
e n E 
e o 

2m e v t 


I kn « v t 6 ett (T « 


T ) 

n. 


+ R + 


s.e . 
i i 


The use of this expression to determine 6^ has two major disadvantages: 
(1) the dc electric field cannot penetrate the plasma without a significant 
associated current, which results in local heating or breaking -down of the 
gas and changes the homogeneity of its properties and (2) the energy-loss 
factor is coupled with the square of the total collision frequency and it 
becomes difficult to separate these two transport parameters and prevent 
uncertainties in v^. from, affecting the determination of & e ££ ■ Both of 
these difficulties can be avoided by using an rf electric field. 

By using equation (16) and the methods of Ginzburg and Gurevich 
(ref. 7, pp. 119 ff. ) it is now poosible to show that when B = 0 in a high- 
frequency electric field (E = E q cos wt , co h 2 irf ) the electron energy 
balance reduces to 


_ _ e 2 n v. E 2 (l +e)R - . „ 

J . E = e ; — r- = 2 kI v,WV T „> + R + 1 Vi (16b > 

2m u + vj i 

e t 

Therefore the difference between the electron temperature T g and the 
reference gas temperature T^ in the absence of radiation (R = 0) and 
secondary eneigy exchange processes (2 s^c . = 0) is given by 


T - T 
e n 


e 2 E 2 (1 + e) n 
o' ' 

3k6 eff m e (tjZ + v t } 


(17) 


When R and 2 s.e. are not negligible we obtain 
i 11 


T - T 
e n 


e 2 E 2 (1 +e)n , . 

, r ~r ■ suAg -r lR + l Vi 1 

3k6effm e (co + v t ) e t eff l 


(17a) 


Generally, when T g - T n * 0 in the absence of an electric field, a 
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steady-state temperature difference is sustained by radiative energy losses 
or secondary energy exchange processes. Equation (17) is valid only if the 
effects of the electric field dominate. 

The quantity e (which arises from the time -dependence of the 
electric field and is given explicitly in reference 7) is usually much less 
than unity when cj » . The quantity II arises from considering the 

frequency and electron-temperature dependent higher order approximations 
to the kinetic coefficients and is also given explicitly in reference 7. It 
is always very close to unity when (a) the predominant collisions are elec- 
tron-neutral particle collisions or (b) the predominant collisions are 

Z 2 

electron- ion collisions but at the same time u » v . The factor II 
differs significantly from unity (and therefore influences substantially 
the electron temperature) only in the case of collisions with ions at low 
frequencies of the applied electric field (to 5 v ) and also in the region of 
gyro resonance if, contrary to the situation treated here, the magnetic 
field is non-zero. 


Equation (17) immediately suggests a method for the determination 

of 6 or 6 , * If a gas or plasma in which R = S s.e. = 0 is illuminated 
k eff & r iix 

by a high-frequency electric field such that >> its energy-loss 

factor can be obtained independently of the electron collision frequency 

and all other transport coefficients, by measuring the rate of change of 

electron temperature with the amplitude of the illuminating electric field 

while the gas temperature T and also of course the reference or 

k, s 

normalization temperature T , remain constant. (The condition of 
constant T implies that the illuminating electric field must heat the 

K , S 

electrons but not the gas. This is virtually assured if the degree of 

2 2 

ionization is small, especially in a flow system. ) Assuming cj » 
and differentiating equation (17) with respect to the electron temperature 
we obtain 


eff 


2e 2 E 


3km u> 
e 


8E 

o 

!FT~ 


E 

o 

T 


aUn&eff) 

ST 


(18) 
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If P is the power fed into the appropriately constructed and 
terminated illumination cavity (propagating primarily the TE^q mode 
as described in Appendix A), Z is the impedance of the gap (377 £1 for 
vacuum, but see Appendix A for values applicable in our case), and A^ 
is the appropriate cross-sectional area of the beam or gap (in our case 
the area perpendicular to the wave propagation axis as described in 
Appendix A), equation (18) can be transformed to 
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eff 


3km go 
e 


(£) (£■) 


8T e 9«n6 e£f )i 

91 nP ' 5T | 
e 


(19) 


When the electron energy-loss factor remains constant or very nearly 
constant with (as it does, for example, in diatomic gases over narrow 

ranges of where the electron energy remains below the threshold for 

exciting the next higher energy-absorption mechanism in the heavy particle) 
we can write 


6 


eff 



3km go 
e 


Z 




( 20 ) 


A convenient determination of the electron energy-loss factor by 
means of equations (18), (19) or (20) can be carried out by measuring the 

electron temperature T in the cavity as a function of the incident power 

e - 1 
of illumination P and obtaining dT^/dE^ = (dE Q /dT e ) or 

dT e /dP = (dP/dT^)"* . The measurement of P is trivial; the measure- 
ment of T is not. 
e 

The alternative method of determining 6^^ from equation (17) 

by measuring the rate of change of electron temperature with the frequency 
2 2 

go at co » Vj. was also considered. Since, however, it is far more 
convenient to change P by one or more orders of magnitude (to correspond 
to desirable conditions for different gases with 6 ££*s varying by as 
much as three orders of magnitude) than go , this alternative approach 
was for the time being abandoned. 
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In both methods, however, the energy-loss factor is inferred 
directly from measurements of macroscopic quantities. Independent 
determination of such parameters as cross-sections, electron mobility 
and particle densities are not required. 


Physical Constraints Placed on the Experiment 

Choice of a method for measuring must be governed by 

several criteria: it must be compatible with the assumptions made in 
deriving equations (18), (19) or (20) and with the environment, geometry, 
other constraints of the experiment; it must introduce the minimum 
number of additional conditions; it must be simple, reliable and economical; 
it must not interfere with the other measurements nor should the other 
measurements interfere with it; and it must yield the maximum possible 
useful information. However, whatever the criteria for choosing a method 
for measuring , they are overshadowed by the necessity to establish 
first the existence of an electron temperature, i. e. , to ensure a Maxwellian 
thermal energy distribution for the electrons in the gas under consideration. 

We discern therefore a certain hierarchy of conditions and con- 
straints that are essential to the success of the experiment. They arise 
from the basic requirement for the existence of electrons with a well- 
defined temperature in a hot gas and the operational requirements of 
heating these electrons with an rf electric field and measuring their 
temperature. 

The first of these constraints that have to be met is the one that 
assures a Maxwellian electron energy distribution. The criterion for 
Maxwellian electrons is that the collision frequency for energy transfer 
between electrons, should be larger than the total collision frequency for 
energy transfer between electrons and neutrals (ref. 7). The energy-loss 
factor for electron-electron collisions is 6^ = . (This value for the 

energy-loss factor is obtained from classical collision theory for two 
particles of equal mass. The general formula is 6 = 2m m / (m + m ) 

which reduces to 2m /m for m » m and to (1/2) for m = m . ) 

6 K ^ C c K 
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Therefore, the effective collision frequency for energy transfer among 

electrons is t * = i- t * where t * is the effective collision fre- 
e,e,e c e,e e,e 

quency for momentum transfer among electrons. From equation (10) the 
total effective collision frequency for energy transfer between electrons and 
all heavy particles is • Therefore, the ratio of these two 

quantities is given by 


R 


M 


1 -1 
7 %,e 

v t 6 eff 


1 
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-1 

T 

e, e 



-1 

K T e, k 


-1 


( 21 ) 


where equation (7) is used for v. 5 . 

t 

If we use the collision cross-section for momentum transfer between 
charged species given as equation (24) in reference 3, we obtain (k + e) 


R 


M 


6 X 10 7 



( 22 ) 


in mks units. The electron energy distribution is Maxwellian when the 
magnitude of the parameter R^. satisfies the inequality R^. » 5 
(ref. 7). It must be noted that at high electron velocities, i. e, , at the 
tail of the distribution function, deviations from Maxwellian are significant 
even at R^ ~ 50 . On the other hand, at low electron velocities (i. e. , at 
low electric fields) the constraint R^ » 5 can be relaxed considerably 
and in our case we can expect a sensibly Maxwellian distribution even at 


Clearly, the most important result of the constraint imposed by 
R^ >> 5 is that the experiment must be operated in an environment 
where n /n is not less than approximately 10 ^ for monatomic gases 

^6 K 

and 10“ or higher for polyatomic gases. It is worth noting that any 
enhancement of the degree of ionization that occurs because of the eleva- 
tion of T makes it easier to meet this constraint. 


e 

When n^ is specified by some other constraint or design 
consideration, this constraint serves to specify the maximum permissible 
neutral particle density n^ . 

The second of these constraints involves the limitation on the 
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electron number density imposed by the requirement that the illuminating 
high-frequency electric field must penetrate the plasma. The condition 
w > Wp = 2irfp , where is the critical plasma frequency in radians/sec 
given by f = 9 (n )*^ 2 in Hz, when n is in electrons/m 8 , ensures that 

p G 6 

the high-frequency electric field will penetrate the plasma and imposes an 
upper limit on the electron concentration given by 

n g < f 2 /81 = o> 2 /324 tt 2 (23) 

where the frequency of the applied electric field is f in Hz or go in 
rad/sec. Any enhancement of the electron number density that occurs 
because of the elevation of the electron temperature makes it tougher to 
meet this constraint and great care should be taken to make sure it is 
observed, (See also the discussion of the eighth constraint, ) 

The third of these constraints arises from the requirement that 
2 2 

go » v, and imposes an upper limit on the neutral particle density n 
t a 

If we consider the case where only one neutral species is present and only 
electron-neutral collisions for momentum transfer are important, the 
neutral particle density must satisfy the inequality 


2 


n 


a 





1. 46 X 10 -8 


2 

GO 



T 

e 


(24) 


obtained by combining special forms of equations (6) and (12), In equation 

(24) the quantity Q is the collision cross-section for momentum transfer 

e a 

between electrons and neutral particles of species a , is the mass 

of the electron and the mean approach speed between electrons and neutral 
particles (g ea ) Q in equation (12) has been taken to be (g ea ) Q - 
(8k T e /irm e ) 1 / 2 . 

Note that although the first constraint can be interpreted to 
specify the minimum permissible degree of ionization, the second and 
third constraints cannot be combined and interpreted to specify the 
maximum permissible degree of ionization. In fact, the second and third 
constraints must be considered independently and to ensure that they as 
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well as the first constraint are met, the number densities of electrons 
and neutrals must be controlled, as far as possible, independently of each 
othe r . 

From this point, in order to make further progress in the design 
of the experiment we must specify the method for the production of the hot 
gas or plasma with independently controlled electron and neutral number 
densities and we must specify the method of measurement of *he electron 
temperature. 

Since a convenient source of hot gas or plasma of independently 
controlled enthalpy, electron number density and neutral number density 
is an arcjet expanding into a vacuum chamber, it was decided to combine 
such a method of production of electrons and hot gas with the most con- 
venient and widely-used method of obtaining a spatially- re solved electron 
temperature in plasmas of this type — the Langmuir probe. A considerable 
body of experience has accumulated over the years in the use of Langmuir 
probes to measure electron temperatures in plasma flows produced by 
arcjets (e. g. , refs. 12, 13 and 14). A cylindrical probe appears especially 
suitable to the conditions expected in this experiment (ref. 12). The use 
of Langmuir probes for accurate determination of electron temperatures 
is not plagued with the difficulties associated with the use of these probes 
for precise measurement of the electron number densities. Nevertheless, 
use of the data obtained by these probes to establish order-of-magnitude 
values of the electron number densities is a welcome addition to the diag- 
nostic information required to establish confidence in these experiments. 
Finally, the data obtained in each determination of T g by these probes 
can establish at a glance the existence or non-existence of a Maxwellian 
electron energy distribution and thereby provide a check on the validity 
of the measurements and tangible proof that the first constraint is 
adequately met. 

Use of this particular method for obtaining introduces 

additional constraints. Observance of these constraints ensures that the 
data obtained by the probe can be interpreted unambiguously to yield the 
electron temperature. 
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Thus the fourth constraint requires that the Langmuir probe 
diameter must not be larger than the mean free path of the 

electrons. This condition implies that the neutral particle density must 
satisfy the inequality 


n 


a 


1 / (D Q ) 
v p ea' 


and for D = 0. 001 m we must ensure 
P 


n 

a 


10 3 Q 


ea 


( 25 ) 


The fifth constraint requires that the diameter of the probe should 
be much larger than the electron Debye length \ D< This condition implies 

e kT rp 1/ 2 

D , » x d ■ = «■<>* Gr-) < 26 > 

n e e 

e 

in mks units and for D = 0. 001 m we must ensure that 

P 

n » 4. 77 X 10 9 T 
e e 

The constraints prescribed by these two inequalities are not 
particularly hard to satisfy and conclude the demands made on the 
experimental environment by the cylindrical Langmuir probe when used 
to determine T 

e 

All other operational constraints are imposed by the technique 
we have chosen for heating the electrons. 

Thus the sixth constraint requires that the characteristic flow 
time of the plasma across the waveguide must be much larger than the 
excitation or relaxation time for the electron temperature (or, in other 
words, that the time -of- flight of the electrons must be large compared to 
the electron temperature rise time). Observance of this constraint 
ensures that the electrons will indeed be heated in the waveguide since 
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their time-of-residence in it will be much longer than the time required 
to heat them. Since the electron energy excitation time is at least as fast 
as the relaxation time 1 / e f£ (ref. 7), if the flow velocity of the 
plasmajet is U and the length of the waveguide in the flow direction is a , 
this condition is satisfied when 


a 

TT 


> 


1 

V. 6 rr 

t eff 


(27) 


This constraint imposes either a lower limit on n (or a minimum per- 

a 

missible length, a , or a maximum permissible flow velocity U ) given 
by 


n 

a 


> 


U 

a ^eff^ea ( ^ea ^ o 


(28) 


where (g ) = (8kT /urn )^^ = 6. 213 X 10^ T . When this con- 

straint is satisfied, we can be certain that at the probe the electron tem- 
perature corresponds to the electric field E q and equations (18) or (19) 
apply. Since the dimension a and the flow velocity U are specified by 
waveguide and arcjet design considerations, this inequality also sets a 

lower limit on n , which is usually compatible with the requirements 
a 

placed on n^ by the first and third constraints. 

The seventh constraint arises from the need to obtain an elevation 

AT g = T - T e j of the electron temperature above some inital electron 

temperature T^ ^ that is equal to or larger than a certain minimum 

(AT ) . . This value of (AT ) . is usually imposed by considerations 

that have to do with the precision of the measurement of T^ . Clearly, 

for example, we must obtain a AT g that is larger than the minimum 

detectable electron temperature elevation with the Langmuir probe 

apparatus used in the experiment. Moreover, this AT^ must be higher 

than the probable error in the measurement of the electron temperature to 

ensure accurate results (see Appendix B). This constraint on a minimum 

allowable AT imposes a lower limit E . on the amplitude E of 
e mm o 

the electric field. If we neglect the dependence of 6^ on T e , we 
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obtain from equation (17) for » v^. an expression for E . in terms 

of (AT ) . given by 

N e mm & 7 


E . 
min 


(■ 




1/2 


T . + (AT ) . - T 

el e mm n 


1/2 


= 3. 83 X 


10" 8 6 


eff 
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el 


(AT ) . - T 

' e 'mm n 


1/2 


(29) 


The equivalent power increment (AP^.^ required to produce (AT ) ^ 

is given by equation (19) as 


1 R ? 

(AP) . = 1. 47 X 10‘ 13 w 6 ,, {-rL) (AT ) . 

mm eff 4Z v e rm 


min 


-4 10 

For example, when 6=10 , w = 1. 54 X 10 rad/ sec = 2, 45 GHz and 

(AT ) . = 100 °K we obtain E . = 58. 9 V/m. Observance of the 
e mm mm 

seventh constraint is ensured when 


E = 3. 83 X 10" 8 6 2 w (T , + AT - T ) 1//2 S E . 

o eff el e n mm 


(30) 


We repeat here that the precise value of the required (AT^)^.^ and 
therefore also of depends on the accuracy in measuring T^ 

obtainable with the apparatus and on the precision required in the deter- 
mination of (see Appendix B). 

The eighth constraint requires that the effect of the high frequency 

electric field on the plasma should be kept sufficiently small so that 

equation (17) remains valid and/or the influence of the electric field on the 

experiment (including the measurement of T e ) remains easily calculable 

(refs. 1 and 7). For example, if the rates of ionization are such that, 

while the gas is in the test section, the production of ion pairs is negligible, 

neither the validity of equation (17) nor the second constraint are violated 

for almost any practical value of E q at f ^ 2. 45 GHz. However, suppose 

E increases above the value of the characteristic plasma field E given 
o P 

in the first approximation by 
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(31) 


Then it is possible, depending on how well the conditions of the first con- 
straint are satisfied, that the electron velocity distribution function deviates 
from Maxwellian, especially at the high-energy tail (ref. 7). To avoid any 
such possible complications it is sometimes advisable to operate the 
experiment so that the amplitude of the applied electric field satisfies the 
condition 


E < E (32) 

op ' 

This constraint is not rigid for all conditions and should be used only when 

-4 

there is danger of violating the first constraint. For example, when . 6 eff = 10 

and co = 1.54X 10* rad/sec, we obtain E =5. 9 T / V/m and co > v? . 

p e t 

Therefore, an electric field amplitude equal to or less than a few hundred 
volts per meter is indicated in cases where there is any difficulty in 
establishing and/or maintaining a Maxwellian electron velocity distribution. 

2 2 

A combination of equations (30), (3i) and (32) gives for u> > 


T . + (AT) . - T 


L “ 


e mm 
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1/2 E q 

" h m ~ (T el + AT e ) 

3. 83 X 10«6 eff 1/Z 1 


1/2 


(33) 


and it is clear that although the seventh and eighth constraints can always 
be made compatible by proper choice of (AT ) . and E , an attempt 

c mill u 

must be made to carry out the experiment with as small an E q or a 

AT /T _ _ (or, at least, as small a AT ) as possible. The additional 

e e, iii — u e 

advantage of this condition, i, e. of operation as close as possible to 

(AT ) . , is that it allows us to observe the variation of 6 ,, with T 

' e min ’ eff e 

in fine detail. 
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When the experiment is operated under conditions that satisfy 
the eighth constraint (i. e. , with electric fields less than a few hundred 
V/m) one would expect that the additional constraints of operation at an 
electric field that is smaller than the breakdown potential of the gases 
involved (ninth constraint) and smaller than the critical electric field 
above which binary theory is questionable (tenth constraint), are also 
satisfied (refs. 1, 7 and 15). 


The eleventh constraint arises from the requirement that the 

elevation of the electron temperature by the applied electric field should 

not be so small that any energy losses or gains experienced by the electron 

due to the existence of small but finite terms of the form 9R/9T and 

e 

9Ss^e./9T e are significant compared to the term (3/2)kn g Vj6 e ^ , while, 


at the same time, AT g should not be so high that the variation of 


eff 


with T is masked. If we define 9 $/9T i-9ft/9T +92s.e./9T we 
e 

observe from equation (16b) that, for to » v* , 
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/ e t o\ o 
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and therefore the first part of this constraint yields 


AT » 
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3 k n v. 6 et 
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(35) 


This inequality implies that the major part of the ohmic heating goes into 
maintaining the electron temperature elevation AT g . 

The second part of this constraint requires that the fractional change 
of 6^ due to the rise of electron temperature should be small f see also 
remarks following equation (19)]. Thus 
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and therefore the eleventh constraint imposes a desirable range of AT g 
given by 


2 A ft 

3 k n v .6 ££ 
e t eff 


« AT 
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r»<* n6 e|f> 

ST" — 


( 37 ) 


The first and second parts of this constraint are compatible except possibly 
in the rare cases when A 9 and 3(ln & e ff)/^T e are both very large. Such 
cases have yet to be encountered. The first part of this constraint can be 
incompatible with the seventh constraint and care should be taken if such 
is the case, to allow for the influence of the terms $ and d9/dT in 

V 

equations (17) and (18). It should be remembered that with the convention 

that s^ is positive, the quantity e. is negative when the electrons gain 

energy and positive when they lose energy by the processes inherent in 

2 s.e. . Thus when, for example, the electrons gain energy by some 
i 

process s^ the values of 8 ^ calculated from equations (18) to ( 20 ) 
without corrections for the influence of 3$/3T e , may be too low. 

Observance of these eleven constraints as well as of the conditions 
leading to equations ( 16 ) to ( 20 ) is necessary for the design of successful 
experiments to determine 6 fi ££ by this method. It is vary fortunate indeed 
that these constraints and conditions are mutually compatible. Conversely, 
these constraints and conditions effectively specify the design of the 
experiment. 


There are two additional important relationships that are essential 
for the design of this experiment. They consist of the expression for 
ohmic heating £2 of a plasma in the absence of a magnetic field and the 
appropriate relationship between the electric field E^ and the illumination 
power P (see Appendix A). 


The first of these expressions is given by 
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in watts per m , and the second by 

E = (4 P Z A )* 2 (39) 

o r 

Under the conditions of the experiment (see Appendix A) the impedance Z 

is close to Z = 420 Q and A = 0. 1016 X 0. 1524 = 1. 548 X 10' 2 m 2 

r 

and therefore 

E = (4Z/A ) 1/2 P 1/2 = 330 P 1/2 (40) 

o r 

where P is in watts and E in V/m. 

o 1 

The operating constraints are summarized in Table I. We can now 
proceed to define an operating envelope for the experiment. 

Operating Conditions 

Practical requirements such as the availability and cost of rf 
power generators, the accessibility of the region within the waveguide, 
the separation between the waveguide walls that is necessary to admit a 
plasmajet of several centimeters diameter and the transmission losses 
associated with ultra high frequencies, resulted in a choice of 2. 45 GHz 
cw for the frequency of the rf power generator. The actual operating 
envelope for the experiment can now be delineated by making use of the 
theoretical considerations and constraints, equations (21) to (40), of the 
previous section. The operating conditions are summarized in Table I. 

As soon as the frequency f is specified, we observe from 
equation (23) that to satisfy the second constraint we must operate the 
experiment with n^ < 7. 42 X 10^ elections per m^. To satisfy the 
first constraint [equation (22), R xf ^ 25 ] for nitrogen with an electron 
number density n g = 7 X 10 and an electron temperature in the jet 

(independent of n ) of T = 3600 °K (at which temrerature Q ~ 

20 e e 4 22 e - 2 
9X10" and 6 „ =5X10 ) we require n . 1 . 2 X 10 nitrogen 

molecules per m . Thus if the static temperature in the jet is 300 °K, we 


25 


require a static pressure in the free jet (approximately equal to the 

2 

pressure in the vacuum tank) of 500 microns (50 newtons/m ) or less. 

To satisfy the first constraint (R.^ — 25 ) for argon with an electron 

number density n = 7 X 10 ^ and an electron temperature o* 3600 °K 

e -20 -5 

(at which temperature Q .-1.2X10 and 6 = 2. 7 X 10 ) we 

24 e > ^ 

require n <1.6X10 argon atoms per m . Thus if the static tempera- 
a 

ture in the free et is 100 K we require a tank pressure of 22, 000 microns 
or less. Therefore, for strict observance of the first constraint the 
vacuum tank should be operated at pressures in the range of 100 to 500 
microns. However, since we have a quick and easy way of determining 
the seriousness of any possible deviations from Maxwellian distribution 
(by observing the linearity of the semilogarithmic plot of the Langmuir 
probe current- voltage characteristics at low currents) we can probably 
extend this range up to 1000 microns. Nevertheless, it should be much 
easier t ^ obtain Maxwellian distributions with argon than with nitrogen 
(and therefore the semilogarithmic plots of the uangmuir probe character- 
istics should remain linear over a wider range of bias voltage for argon 
than for nitrogen). 

23 

The third constraint, equation (24), requires that n^ < 3. 4 X 10 

n 24 a 

for nitrogen at T = 3600 U K and n < 2. 5 X 10 for argon. These 
e a 

numbers are compatible with the results of our effort to comply with the 
first constraint, i. e. , if we satisfy the first constraint we also satisfy 
the third. 

When a Langmuir probe of 1 mm diameter is used, the fourth con- 
straint, equation (25), requires that the neutral particle density n be 

22 

equal to or less than 10 (corresponding to 400 microns tank pressure at 

2 3 

a static temperature of 300 °K ) for nitrogen and 0. 8 X 10 (corresponding 
to 1100 microns tank pressure at a static temperature of 100 °K ) for 
argon. This constraint also is compatible with the first constraint. 

13 

The fifth constraint, equation (26), requires that n^ » 5 X 10 
and this is compatible wrtfi the second constraint, i. e. , we can satisfy 
the second constraint while we also satisfy this constraint. 


26 



With a = 0. 1524 m and U ^ 3000 m/sec, the sixth constraint, 

2 1 

equation (27), requires that n >1.2X10 for nitrogen (with T = 

2 1 a e 
3600 °K) and n > 9 X 10 1 for argon (at the same T ). This constraint 
a 6 

is also compatible with the first constraint and implies that the pressure 
in the vacuum tank should remain above approximately 40 to 100 microns. 

The seventh constraint, equation (30), requires that the applied 

electric field E q be greater than the minimum electric field required to 

obtain a minimum elevation of the electron temperature (AT e ) m ^ . Now 

assuming that can be determined within ± 100 °K and that we want 

the uncertainty in to be of the order of 10% , we conclude that the 

minimum required electron temperature elevation is (AT ) £ n = 1000 °K . 

To obtain this elevation at T 4 = T requires an electric field of approxi- 

e 1 n 

mately 418 V/m (4. 18 V/cm) or 1. 62 watts of illumination power with 
and 99 V/m (0. 99 V/cm) or 0. 090 watts with argon. 

Assuming the initial electron temperature T^ _ q 2600 °K 
while AT^ = 1000 °K, we obtain E^ = 795 V/m for nitrogen and 188 V/m 
for argon. The eighth constraint, equation (32), requires that E q ^ E^ 
and we see that it can be satisfied, if necessary, without much difficulty. 

The ninth and tenth constraints are also satisfied when the eighth 
constraint is satisfied. 

A careful estimate of the magnitude of the term 3$/ OT^ from 
the data provided by Allen (ref. 16) and Morris et al. (ref. 17) indicates 
that the first part of the eleventh constraint, equation (35), is satisfied 
when AT g > 50 °K for nitrogen and AT g > 1000 to 2000 °K for argon 
(depending on the enthalpy). Note that the presence of large concentrations 
of metastables in the argon jet at certain plasmajet enthalpies may require 
the operation of the experiment at AT considerably higher than this 
lower limit. On the other hand, if the concentration of meta stable argon 
in the jet is somehow suppressed, the required AT^ may be considerably 
lower than this limit. For practical purposes the effect of 3§>/3T e should 
be investigated in argon experiments by observing the variations of & e f£ » 
as determined by equations (19) or (20), with plasmajet enthalpy and 
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applied AT^ (or illumination power) and care should be taken to carry out 
the measurements in a region where the variation of & e ££ with these two 
parameters is negligible. 

An estimate of the magnitude of the term 9(in ^ e £f)/^^ e ^ rom the 

room-temperature data of Crompton and Sutton for nitrogen (ref. 18) shows 

that even at ~ 11, 000 °K (where begins to vary rapidly with ) 

the magnitude of this derivative is such that [ 9(i n ^ e ££)/^T e ] * « 

&eff AT e /AS e ££ ~ 2200 °K and therefore an elevation of the electron 

temperature by any amount less than 2200 °K will satisfy the second part 

of the eleventh constraint, equation (36), since provided we operate below 

T = 11, 000 °K (say, at less than 10, 000 °K) any elevation by an amount 

less than 2200 °K will automatically be much less than dT^/d(& n 8 e ££) 

evaluated at T < 11, 000 °K (AT /A5 decreases rapidly at T ^ 11, 000). 

Therefore, for nitrogen we can write for the required electron temperature 

elevation, A T 

9 e, req 

50 °K < AT < 2200 °K, for N 0 , T < 10,000 °K 

e, req ’ 2 e ’ 

where the "much smaller than" requirement of equation (36) has been 
relaxed since the term 6 e ^AT e /A6 e ££ is considerably larger than 
2200 °K at T e < 10, 000 °K. For argon, by definition, 5^ = 6^ is 
constant with T^ and therefore the second part of the eleventh constraint 
is always satisfied. The electron temperature elevation required for 
determination of fine structure in 6 ££ = ^ e ff (T e ) f° r argon is therefore 
primarily dependent on the magnitude of the masking effects due to 
3$/9T e . These effects can be kept small when 

AT > 1000 to 2000 °K, for A, T < 7, 000 °K 

e, req ’ e 

It should be noted that, for cj = 1. 54 X 10^ rad/ sec, the ohmic 
heating term given by equation (38) becomes, with the aid of (40), 

6 = 5.93X10‘ 29 n vT 2 = 6. 42X10" 24 n v, P (41) 

e t o e t 
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At = 3600 °K and the conditions imposed by the first and second con- 
straints, we obtain 

n = 1.81X10‘ 4 P, W/cm 3 , for N 2 

and (using typical momentum cross-sections for argon given by reference 
19) 

it = 3.22X10 ' 3 r, w/cm 3 , for A. 

We observe from these values of fl and the power requirements 
of the seventh and eighth constraints that the increase in gas enthalpy due 
to ohmic heating of the electrons will indeed be small even if we assume 
instantaneous electron energy relaxation. 


Experimental Apparatus 

Vacuum tank. - The experiments are carried out ii. a 4-foot 
(1. 22 m) diameter, 10-foot (3. 05 m) long stainless- steel vacuum tank with 
the arcjet mounted on the 4-foot diameter access door that also serves as 
the front bulkhead or closure. The axis of the arcjet is at the center of 
the door and coincides with the axis of the tank. The primary vacuum 
system consists of a Rootes-type dry blower (Stokes Model 1722-S Mechani- 
cal Booster Type high- vacuum purnp with 1600 cfm maximum displacement 

-4 

and an ultimate blank- off pressure of below 2X10 torr, driven by a 
10 hp open motor) backed-up by a Stokes Microvac Model 412-H water- 
cooled rotary-piston mechanical high-vacuum fore pump (with 300 cfm 
displacement, driven by another 10 hp motor). Both of these pumps were 
manufactured by the Stokes Equipment Division of Pennsalt Chemical 
Corporation. With this vacuum system the tank can reach a pressure of 
1 micron in less than 20 minutes of pumping time. The tank is equipped 
with an 8 -inch high- vacuum gate valve that can be used to throttle the 
vacuum pumps and a 1-inch ball valve that can be used to bleed air into 
the tank. The tank is connected to the vacuum system through an 8-inch 
diameter flexible bellows. 
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The tank is equipped with five 12 -inch portholes and one 6- inch 
porthole that serve either as viewports or provide convenient locations 
for feedthroughs for the various gas-, power-, and water-lines that are 
connected to apparatus inside the tank. The space between the front and 
rear 3 -foot-high legs of the tank is filled with a one -foot deep slab of 
concrete. This slab weighs in excess of 4000 lbs and was poured around 
several large bolts set into the concrete floor. It has proved effective in 
reducing the transmission of vibrations from the vacuum system to the 
tank and the equipment in it. 

Arcjet system. - The arcjet system consists of a Thermal 
Dynamics Corporation 50N torch provided with a 0. 65 -inch long, 0. 625- 
inch diameter plenum chamber downstream of the arc. The plenum 
chamber exhausts into the vacuum tank through a conical supersonic 
nozzle of 8° 30 1 half -angle with a throat diameter of 0. 165 inches and an 
exit diameter of 0. 500 inches. The diameter of the hole in the anode 
insert used with this nozzle is 0. 188 inches. For lower flow rates (less 
than 0. 5 gm/sec) a supersonic nozzle of the same half-angle but a throat 
diameter of 0, 091 inches and an exit diameter of 0. 425 inches can be used. 
The anode insert diameter for this nozzle is 0. 100 inches. This arcjet is 
powered by three TDB1A-14 selenium rectifiers manufactured by the 
Miller Electric Manufacturing Co. Each of these rectifiers is rated at 
14 kW on a continuous duty cycle and 20 kW maximum output. 

The primary gas is injected into the arc chamber of the arcjet. 

A secondary gas can be injected into the plenum chamber. The primary 
and secondary gases are metered through calibrated tape red -glass -tube 
meters. Up to four different separately-metered gases can be mixed in 
any desirable proportions before injected into the primary or secondary 
gas feed lines. Four Simet Proportional Mixers (each consisting of a 
throttling valve and a flow meter and made by National Welding Co. ) are 
used for the injection of these additive gases. Besides the primary and 
secondary flow meters the arcjet control console is equipped with a volt- 
meter, an ammeter, primary and secondary gas flowmeter pressure 
gauges, hf arcjet starter, off- and- on control, arcjet power control, water 
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coolant flow meter and water coolant entry and exit thermometers. The 
primary and secondary gas temperatures are measured by additional 
thermometers placed upstream of the gas flow meters. 

This arcjet can easily attain enthalpies of up to 40 MJ/kg 

(17, 200 BTU/lb) and stagnation temperatures of approximately 8, 000 °K 

with nitrogen at plenum chamber pressures (stagnation pressures) in the 

vicinity of 1 atmosphere. It can attain considerably higher stagnation 

temperatures with argon or helium (approximately 14, 000 °K). The 

manufacturer (Thermal Dynamics Corporation) claims that it is 

possible to attain enthalpies up to 260 MJ/kg (60, 000 BTU/lb) with 

nitrogen (with corresponding stagnation temperatures between 10, 000 

and 14,500 °K indicating that the stagnation pressure ranges between 
- 4 - 1 

4X10 and 2X10 atmospheres) and stagnation temperatures up to 
28, 000 °K with argon and helium. Under the usual operating conditions 
the arcjet chamber pressure is in the range of 0. 1 to 1 atm with nitrogen 
and 0. 5 to 2. 5 atm with argon. 

Microwave illumination system. - A schematic of the microwave 
system used to make the necessary measurements is shown in figure 1. 

The Raytheon Corporation rf power generator is stable when operating 
above 10 watts cw and for this reason, rather than operate the unit at very 
low power levels, an attenuator consisting of a long coaxial cable in series 
with a variable attenuator is used to give stable operation in the lower 
ranges of illumination power. Figure 2 is a schematic of the special wave- 
guide, designed to illuminate the plasma on the basis of the calculations 
given in Appendix A. The electric field in this waveguide test section has 
been measured with a carefully calibrated search coil (traversed down the 
centerline of the waveguide, i. e. , the x-axis) and both its magnitude and 
distribution have been found to correspond almost exactly to the theoreti- 
cally predicted field that one expects for a TE^q dominant mode (see 
Appendix A). Since the experimentally determined E q was smaller by 
15% or less from the theoretically predicted value of E q while all the 
assumptions made in reducing the experimental data on the determination 
of E q were such that they are expected to lead to slightly smaller values 
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of E q , the theoretical value of the electric field given by equations (39) 
and (40), was used in interpreting all experiments. 

The microwave illumination system is equipped for measuring 
(a) the power fed into the waveguide, (b) the power leaving the waveguide 
and (c) the power reflected from the waveguide. Since the waveguide and 
associated connector losses were found to be symmetric, i, e. , they were 
the same no matter from which end of the waveguide the power entered, 
the illumination power P is determined by summing the input and output 
power and dividing by two. Usually the illumination power is 80% to 90% 
of the power input to the waveguide. 

Measurements of the reflected power are carried out to ensure 
that reflection is negligible (and therefore the second constraint is observed) 
and that the insertion of the Langmuir probe produces negligible or no 
interactions with the rf field in the waveguide. 

Langmuir probe system and traversing mechanism. - Electron 
temperature measurements are carried out with 0. 2 and 1 -mm-diameter 
cylindrical Langmuir probes. These consist of a tungsten wire of the same 
diameter inserted through a quartz capillary tube with O. D. of approxi- 
mately 0. 5 to 0. 7 cm and I. D. slightly larger than the O. D. of the wire. 

The front end of the quartz tubes is given a conical shape of approximately 
30° half-angle and the tungsten wire is pulled through a hole at the center 
so that it protrudes the desired length ahead of the cone (5 to 7 mm for the 
1 mm diameter wire). At the lower enthalpies (less than 20 MJ/kg) these 
probes are effectively cooled by radiation and so far they have never 
melted, thanks to the low pressures in the vacuum tank that decrease the 
heat transfer rate and the enthalpy flux. The tungsten wire is long enough 
to reach outside the tank through a special rubber-and-vacuum-grease seal. 
Since the distance from the probe tip to the seal is approximately 8 feet, 
the wire at the seal is never significantly above room temperature. 

A gas-cooled probe was also designed and tested for operation at 
higher enthalpies. This probe consists of an outer jacket made of a 
seamless tube of 90% platinum- 10% iridium with an outside diameter of 
0. 050" (1. 25 mm) and a wall thickness of 0, 007", over an inner 
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hypodermic 304-stainless steel tube with an outside diameter of 0. 028" and 
a wall thickness of 0. 004". The outer Pt-Ir jacket is joined to an outer 
0, 187" O. D. 304-stainless steel tube, of 0. 022" wall thickness, through 
a flare adaptor. The flare adaptor is soldered to the stainless outer tube 
and the platinum jacket of the probe with Easy Flow 45 solder. Inside , u he 
0. 187" O. D. outer tube is a 0. 125" O. D. 304-stainless steel tube that is 
soldered to the inner hypodermic 0. 028" O. D. tube of the probe. The 
coolant gas comes in, flows through the 0. 125" tube, squirts through the end 
of the 0. 028" hypodermic against the front closure of the Pt-Ir jacket and 
returns through the space between the hypodermic and the jacket and 
between the 0. 125" and 0. 187" tubes and finally through a hose connected 
to a feedthrough at the vacuum tank wall, to be exhausted outside the tank. 
The gas is supplied, at 400 lb/in , to the 0. 125" O. D. inner tube through 
a high-pressure hose connected by means of another feedthrough to a 
compressed -nitrogen cylinder outside the vacuum tank. All connections 
and feedthroughs are carefully tested for leaks before each experiment. 

The 0. 187" O. D. tube is inserted into a 0. 7 mm O. D, quartz tube whose 
tip is tapered so that it fits closely around the 0. 050" O. D. probe which 
protrudes 7 mm beyond the quartz. 

Although some initially promising results were obtained with this 
probe, it failed, for reasons as yet undetermined, at the point where the 
Easy Flow solder joins the platinum tube to the flare adaptor that connects 
the platinum probe to the outer 0. 187" O. D. stainless tube. 

The Langmuir probe is mounted on a three-degree of freedom 
traversing mechanism. The traversing mechanism rides on two 9-foot- 
long 1-1/4" square, hollow cross-section, stainless- steel rails, water 
cooled, that are mounted firmly on the right and left walls of the tank 
parallel to the tank axis and about 3 feet from the bottom of the tank. 

An outer carriage rides on eight ball-bearing rollers that ride 
backward and forward on the rails to provide traversing in the flow direc- 
tion (x-axis of figure 1). The outer carriage contains two more rails, 
approximately 30" long, mounted in such a way as to provide traversing 
in the right-and-left direction (z-axis of figure 1). An inner carriage 
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rides on these rails and contains a worm-screw gear and two stabilizing 
rails for up-and-down motion of the traversing platform. Mounted on 
this platform is a round heat-exchanger and coaxial with the heat- 
exchanger and through a small hole at the center of its backplate is 
mounted the 22-inch- long quartz probe shield surrounding the stem of the 
probe. The Langmuir probe usually protrudes about 15 to 18 inches in 
front of the backplate of the heat-exchanger. The end of the quartz probe 
shield is usually 4 to 7 inches behind the heat-exchanger backplate and is 
therefore completely shielded from the hot gases. 

The length of the quartz tubes that form the outer shield of the 
stem of the probe is 22 inches. 

Back-and-forth motion (x-axis) is provided by continuous wire 
belts attached to both sides of the outer carriage parallel to the rails and 
driven by a high- vacuum -rated electric motor mounted inside the vacuum 
tank. Right-and-left motion (z-axis) is provided by a similar drive 
mechanism mounted on the outer carriage. Up-and-down motion (y-axis) 
is provided by a high- vacuum -rated electric motor that drives the worm- 
screw. The position is sensed by Selsyn motors (motion and position 
transmitters and receivers) and Veeder-Root counters (position displays). 
Thus only electrical power lines enter into the vacuum tank to actuate and 
readout the traversing mechanism. The absolute accuracy of this system 
is ±0. 005" and tank warping is thought to be the limiting factor. 

The heat-exchanger is so efficient in removing the energy from the 
gas that neither the tank walls nor the pump inlet experience a temperature 
that rises significantly above room temperature. 

The Langmuir probe and the heat-exchanger are rigidly mounted 
together and move simultaneously and predictably so as to traverse 
approximately 6 feet along the length of the tank (x-axis or flow direction) 
approximately twenty-four inches right-and-left (z-axis) and sixteen inches 
up-and-down (y-axis). The position displays and traverse controls are 
mounted in a control panel outside the tank. 
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The electrical circuit for the Langmuir probe is substantially the 
same as used by Kelly et al. (ref. 12) except that their 100 ft resistor, 
which is used as a shunt to measure the current, is replaced in our 
apparatus by a high-precision (better than +0. 5% of full scale) decade 
capable of varying the va i<» r \his resistance from 1 to 99, 999 ft . The 
maximum bias voltage for the probe is 24 volts (supplied by four dry cell 
batteries, such as Eveready 510S, connected in series, each battery made 
up of four size F cells). A four-position switch is provided that gives zero 
bias (off position), full positive bias (zero to +24 volts), negative and 
positive bias (-12 to + 12 volts) and full negative bias (zero to -24 volts) 
on the probe. The probe bias in each switch position is varied by a coarse 
sweep (provided by a one-turn dual pot with 25, 000 ft in each section rated 
at 2 watts for each section, the two sections series-connected for a total of 
4 watts) and a fine sweep (provided by a one -turn potentiometer with a 
resistance of 5000 ohms and a rating of 4 watts). Ten-turn pots were also 
used for these sweeps but after some experience was obtained they were 
replaced by the simpler pots. As the potentiometer sweeps the bias 
voltage on the probe, it is recorded on the x-axis of an x-y recorder 
(Monels 135 -A or 7000- A, Moseley Division of Hewlett-Packard Corp. ) 
which also records the voltage drop acrss the decade (usually set at 
100 ft ) on the y- axis, to obtain the current-voltage characteristics of the 
Langmuir probe. 

Probe noise spectrum analyzer and integrating readout. - As a 
preliminary step towards developing an automated noise -suppressing 
electron-temperature (ANSET) readout system a Hewlett-Packard 3460 
integrating voltmeter with 0. 1 sec maximum integrating time was used 
to make a preliminary analysis of the noise spectrum of the Langmuir 
probe operating in a nitrogen plasma in order to decide (a) if the noise is 
truly random and can be averaged or integrated out and (b) what is the 
optimum averaging time for improved signal-to-noise ratio. The results 
of more than 100 different readings at each of several probe bias voltages 
were plotted as an integral spectrum and the RMS noise and percentage 
noise were obtained. It was concluded that (a) the major source of noise 
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seems to have Gaussian distribution, (b) an averaging time of 10 seconds 
would ensure that probe currents can be measured to the order of better 
than 1% of their value at any probe bias setting, (c) the RMS noise current 
seems to be roughly proportional to the square root of the electron current, 
and at the highest electron current, where the RMS error is about 4% , one 
can make a model of the electron current coming in bursts and calculate 
that the bursts' frequency should be about 6 kHz. An examination of the 
probe current by a fast oscilloscope indicated that, under some arcjet 
operating conditions, current bursts do indeed appear to be coming at about 
this frequency. 

That the source (arcjet) is stable for 10- second intervals was shown 
by using an x-y plotter to obtain time recordings of the current to the probe 
at several bias voltage settings. During the course of these measurements 
it was discovered that there seems to exist a long-term drift in the arcjet 
current and voltage charcteristics that causes the probe potential to drift 
or change. These drifts show up as a sudden (i. e. , in less than 0. 1 sec) 
drop of the arcjet voltage by 1 to 3 volts followed by a slow but accelerating 
drift of the voltage to approximately its original v' lue. The period of this 
drift is usually of the order of 10 minutes. While this drift occurs the 
probe characteristics also drift (i. e. , are displaced) on the bias voltage 
axis without a change of the shape or slope of the curve. 

The results of these studies indicate that this arcjet can operate in 
two distinct modes: The "noisy arc" mode (where all evidence points to an 
intermittent arc that strikes and is extinguished with a period of the order 
of 100 )jl sec) and the "quiet arc" mode (where the arc operates continuously 
while the plasma and the probe signal are quiet). This quiet arc mode is 
encountered when the arc chamber pressure is of the order of 0. 15 atmos 
absolute (in nitrogen). The probe then is so quiet that the noise is negli- 
gible. Drifts occur with much lower frequency when the arc operates in 
the "quiet" mooe. However, not all operating conditions and requirements 
can be satisfied by arcjet operation in the quiet mode and the effort to 
develop an integrating probe readout that enhances the signal-to-noise 
ratio was continued. 
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As a preliminary step, before employing digital integrators in the 
ANSET readout system, it was decided to extend the linearity of the 
Langmuir probe current-voltage plots by a full decade by surveying very 
accurately the region close to 0 volts (i, e. , to establish the linearity down 
to the 0.1 mv/inch range of the recorder whereas previously it was limited 
to 1 mv/inch). The probe resistance remained at 100 ft. Data were taken, 
based on a less "^'irate integration procedr " as follows: an automatic 
voltage steppe ouilt and was connecte -wo x-y recorders, one of 

which plotted age applied” (in steps) vs “time”, and the other plotted 

the corresponds g ”probe current” (voltage developed across a 100 ohm 
series resistor) vs "time”. By a procedure of visual integration of the 
graphs, the above technique allows a fairly accurate estimate of the expected 
value (mean) of the stochastic process describing the current through the 
probe. Since the power -density spectrum of the noise was found to have 
all components in the range of 10 Hz to 10 kHz, it was concluded that an 
integrating time of 5 sec per voltage -step is long enough to minimize the 
noise effect, and short enough to keep drift from becoming a problem. 

The automatic -voltage stepper* s versatility also permits us to 
concentrate on the most critical region of the graphs and to specify 
accurately the values of the coordinates there, by essentially magnifying 
this region ten-fold. 


Repeated measurements of 5^ with this technique yielded values 
very close to those previously reported, as was to be expected since the 
” integrating” was again done by the x-y plotter and the human eye. Never- 
theless, the linearity of the probe characteristics was extended over an 
additional decade and useful measurements of 6 j*j could be obtained in 
the presence of considerably higher noise levels. 


The automatic voltage steppe^ (AVS), is fast and uses integrated 
circuits, of the type used in the newest digital computers, throughout. 
Timing is controlled by an extremely accurate, temperature controlled 
100 kHz crystal oscillator employing unijunction semiconductors. Front- 
panel controls allow starting and stopping the sequence of steps at any 
time, as well as monitoring the progress of the sequences. The AVS is 
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only one of the four main components of the ANSET readout system. The 
other components consist of a Hewlett-Packard 2401C integrating digital 
voltmeter, the logic circuits coupling the output of the integrating voltmeter 
to a coder (card, paper or tape punch) and the computer and computer 
program that convert the coded data into an electron temperature. 

Testing of this ANSET readout system is incomplete. From the 
available data it appears that the accuracy of this or similar fast and so- 
phisticated digital integration equipment is really necessary before the 
question of the influence of noise on the measurements is given a final and 
authoritative answer. 

The spectroscopic apparatus used to determine the energy state of 
the hot gas by observation of the spectral intensity of the radiation emitted 
from the plasma when it is excited by an electron beam, is described in 
the next section. 


ENERGY STATE OF THE GAS 

To define the energy state of the hot gas of these experiments, an 
electron beam device was built and used to ionize and excite directly rota- 
tional and vibrational bands of diatomic components of the plasma. The 
rotational and vibrational fine structure of the spectrum of the radiation 
emitted from the molecules ionized by the beam was then used to determine 
the rotational and vibrational temperatures of these moledules. 

Sp< ctroscopic analysis of the radiation emitted when a gas is 
excited by an electron beam has been used by several investigators (e. g. , 
Muntz and Marsden, reference 20, Petrie, references 21 and 22, and 
Sebacher, reference 23) to define the energy state of vibrationally and 
rotationally relaxing high-speed flows of nitrogen and other gases. 

The electron beam diagnostic technique is particularly attractive 
because (a) the small beam diameter (typically t mm) and small beam- 
spreading leads to good spatial resolution and requires no symmetry 
(axial or otherwise) of the plasma jet in ^rder to obtain distributions of 
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the observable plasma parameters and (b) it has no measurable influence 
on the basic energy state of the test gas and can be applied to any low- 
density flow regardless of the gasdynamic, chemical and thermodynamic 
nature of the gas. In this technique, described in detail by other authors 
(refs. 20 to 24), the narrow beam of electrons is projected across the flow 
and the interaction of the electrons with gas particles produces a thin 
column of radiating gas that is nearly coincident with the beam. In air or 
nitrogen the emission is excited by collisions between beam electrons and 
nitrogen molecules in the ground electronic state (N^X^S) . In specifying 
the detailed mathematical model for analysis and interpretation of the 
observed radiation, it is assumed that (1) excitation occurs directly from 
the ground state to the N^B 2 state of the nitrogen ion, (2) the 

Franck- Condon principle applies directly in the excitation process and a 
spontaneous radiative transition of the excited ion to its ground electronic 
state (N^ X 2) occurs before any collisional process takes place, 

(3) Born ! s approximation holds (since the energy of the beam electron is 
large compared to the ionization energy of ) and the rotational excita- 
tion process obeys the usual optical selection rules (i. e. , the probabilities 
of transition between the rotational energy levels are assumed equal to the 
corresponding values for the usual optical transitions) and (4) an expression 
for the population in the ^ ground state is available (in this case 
Boltzmann population distributions for the rotational and vibrational tem- 
peratures, references 25 and 26). Previous analyses (refs. 20 to 24) 
substantiate the validity of the assumed excitation process under special 
conditions. Similar conditions that justify these assumptions also prevail 
in our experiments and therefore the rotational and vibrational tempera- 
tures determined by this technique are those of the un-ionized nitrogen 
molecules in the ground electronic state, i. e. , the corresponding tempera- 
tures of the molecules before excitation by the electron beam. 

When an electron beam with energies in the range of 10 to 30 Kev 
is passed through air or molecular nitrogen, the predominant radiation is 
due to the first negative emission system of . The most intense band 

in the emission is the (0, 0) band at 3914 A. Radiation has also been 
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observed from bands in the (0, 2), ^0, 1) and (1, 0) progressions. In 
addition, bands from the second positive system of N ? have been observed 
with intensities near that of the (1,0) band. 

The rotational temperature of is determined by measuring 

L + 

the relative emission intensities of the rotational lines in any given 
band [usually the (0, 0) band at 3914 A] . The bands consist of a 

R -branch and a P-branch of rotational lines. The P-branch lines are very 
close together and form the band head. The R -branch lines are spaced 
rather evenly towards the violet with increasing rotational quantum number 
values and can be resolved by a spectrograph with moderately good resolu- 
tion. The rotational temperature is obtained from the measurement of the 
intensities of the R- branch lines by the straight-line plot method discussed 
in the literature (refs. 20 - 24). The accuracy of this determination of the 
rotational temperature depends on the sensitivity of the optical system used 
(i, e. , on the signal-to-noise ratio) and it appears that it is probably limited 
to ± 1% or less even for a system of moderate sensitivity and resolution. 
Such a system can be expected to achieve a practical accuracy of ±3% 
in the determination of rotational temperatures by these techniques (refs. 

20 -24). 

The vibrational temperature is obtained from the measurement 
of the integrated emission intensities of at least two vibrational bands. 

The bands in both the (0, 2) and (0, 1) progressions of the first negative 
emission system of are suitable for vibrational temperature measure- 

ments. The (0, 1) pi ogression has intensity peaks at 4278 A for the (0, 1) 
band, at 4236 A for (1, 2), at 4199 A for (2, 3) and at 4167 A for the (3, 4) 
band. The (1,0) progression is overlapped by the (0, 1) progression of the 
second positive system. Since the Franck-Condon factors and wave- 
numbers of the bands are known, a measured integrated intensity ratio, 
e. g. , between the (0, 1) and (1, 2) bands or the (0, 1) and (2, 3) bands, can 
be used to obtain the magnitude of the vibrational temperature. The 
accuracy of this vibrational temperature determination depends mainly on 
the technique employed for intensity measurement; however, uncertainties 
in the Franck-Condon factors and errors associated with the integration of 
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the band intensities probably limit the accuracy of this temperature 
determination to approximately ±6%. A system of moderate sensitivity 
and resolution can be expected to achieve a practical accuracy of ±9% 

(refs. 20 - 24). 

Such an optical system, consisting of a Jarrell-Ash 0. 5 -meter 
Ebert scanning spectrometer and an EMI 9502 S photomultiplier tube, was 
built and used to obtain rotational and vibrational temperatures in the 
nitrogen plasmas previously described (ref. 1). The spectrometer used 
a grating with 1180 grooves/mm blazed for 4000 A. The effective aperture 
ratio of the spectrometer was f/8. 6 and the resolution was at least 0.2 A 
in the first order. Typical dark current of the photomultiplier was 0. 5 
nanoamps. The photomultiplier signal was processed by a Keithley 
Instruments Model 417 picoammeter and recorded by a Moseley 7000A 
x-y recorder synchronized with the scanning drive of the spectrometer. 

A separation of approximately 1 cm between the peaks of the Hg 3131 A 
doublet (3131. 55 A and 3131. 83 A ) was easily obtained by the recorder. 

When it was deemed necessary to eliminate or check the drift in the electron 
beam intensity or other sources of noise, zero-shift or poor response, as 
well as the effect of integrated noise, a light beam splitter was used to 
divert part of the light beam to another photomultiplier (RCA 1P28). This 
reference photomultiplier was used to monitor total radiation intensity of 
selected parts of the spectrum (by means of selected light filters) while 
scanning (e. g. , total intensity of the (0, 0) band) and its signal was also 
processed by a Keithley 417 picoammeter. The electrical signal from the 
spectrometer photomultiplier (EMI 9502S) was divided by the signal from 
the reference photomultiplier by an on-line analog computer (Electronic 
Associates TR-20, 1967 model). This division eliminates the effects of 
gradual changes in beam current which might occur during a scan as well 
as other sources of possible error due to time -dependent fluctuations in 
light intensity. The Keithley 417 picoammeter s proved particularly con- 
venient because they include a bucking network so that d. c. noise can be 
eliminated. For calibration of the optical system a standard of spectral 
irradiance traceable to the National Bureau of Standards (coiled-coil 


41 



tungsten filament quartz -iodine lamp calibrated over the spectral range 
from 2500 A to 25, 000 A ) with a special power supply was obtained from 
Eppley Laboratory, Inc. 

The spectrometer light- gathering system had an effective aperture 
ratio of f/4 . The spectrometer slit image was focused on the space occu- 
pied by the electron beam by backlighting. The image was approximately 
1 mm high. Therefore the volume of radiating plasma examined by the 
spectrometer is of the order of 1 mm cubed. The electron beam was at 
right angles to the plasmajet axis and shooting down from the top. The 
spectroscope axis was at right angles to the electron beam axis on the left- 
hand side of the plasmajet axis and almost perpendicular to it, (The angle 
between the optical axis and the perpendicular to the plasmajet axis was 
6.5°). An observation port was on the right-hand side of the plasmajet 
axis. The cylindrical Langmuir probe was on the plasmajet axis (pointing 
upstream) at approximately the same distance from the arejet nozzle exit 
plane as the electron beam (approximately 16 cm). For the electron beam 
experiments the Langmuir probe was moved 15 cm further downstream 
and the microwave apparatus was removed from the vacuum tank. 

For these measurements, the arejet was operated at an energy 
input to the gas of between 2x10^ and 4X10^ joules/kg and the stagnation 
pressure in the plenum chamber varied from 0, 5 to 0. 6 atmospheres, 
corresponding to the conditions of reference 1. 

Table II shows typical vibrational and rotational temperatures, 

T v ib ^rot res P ectivel Y> measured with this apparatus in a pure 

nitrogen plasma. Table II also shows the corresponding free electron 
temperatures T^ ^ _ q zero illumination power and with the 
electron beam off) obtained with the Langmuir probe, the corresponding 
stagnation temperature T obtained from the average enthalpy derived 
from an arejet energy balance, and the corresponding static temperature 
T computed from an isentropic expansion from the stagnation pressure 
and temperature to the vacuum tank pressure of 400 microns. It is seen 
that T g g _ o * n substantial agreement with the measured vibrational 
temperature and ^ rQt is m substantial agreement with T . The 
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implications of these measurements are that under the conditions of the 
experiment the electron temperature is closely coupled and in near 
equilibrium with the vibrational temperature even for pure nitrogen, as 
previously suspected (ref, 1). In addition, the rotational temperature 
appears to be closely coupled to the static (translational) temperature as 
expected for these conditions (refs. 20 - 24). 


Now consider the case of thermal equilibrium between the various 

degrees of freedom of the molecule, under conditions such that the 

radiative loss term R and the secondary energy exchange term Ss.e. are 

i 1 1 

both negligible. Then the contributions to the energy loss factor from the 

various excitation mechanisms (i. e. , elastic or translational, rotational, 

and vibrational) can be simply added to yield the total inelastic energy loss 

factor, 6 or 6 ££ . Moreover it can be shown that under certain 

reasonable assumptions each of these contributions 6. , becomes inde- 

pendent of the associated excitation temperature T7 of the molecule when 

this temperature is much larger than the characteristic temperature 0. 

for that process, i. e. , the contribution 5^ is independent of for 

T.' » 0. and T » 0. . For example, the contributions to 6 from 
J J e j r ' k 

vibrational excitation, 6 , becomes 


b 

v 



1,0 


(T e ) • (1 - e 



1 ) - 


e /(e 

v 


0 /T 
v vib 


1) x 


X (T 
v e 



(42) 


where Q (T e ) is t ^ ie de -excitation probability for the first excited 

vibrational state by electron collisions, 0 v is the characteristic tempera- 
ture for vibrational excitation and = T/ is the vibrational tempera- 


v Equation (42) can be obtained by specializing the theory of translational - 
vibrational relaxation as outlined by Montroll, E, W. ; and Shuler, K. E. : 
J. Chem. Phys. 26, (1957) p. 454, to the case of electrons as collision 
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It is easily seen from equation (42) that when T and 

» 0 v ; is independent of T . Similar arguments can be made 

for the other excitation processes. At the same time 6^ is negligible 

when T « 0. and T. V « 0. . When the gas is not at equilibrium each 
e J ) J * 5 

6. is computed at the temperature TV (ref. 3). Therefore, for the gases 
3 J o 

and in the range of temperatures of interest to us (roughly 300 K to 

10, 000 °K) the contribution to 6 or 6 from rotational excitation 

energy loss mechanisms is substantially independent of the rotational 

temperature and only the contributions from vibrational and electronic 

excitation are dependent on the associated excitation temperatures T 

and T e ^ ec . Similarly for the gases and conditions of interest to us, 

R and Ss.c. can be considered negligible. 

This is a consequence of the fact that 9 rot is 2. 87 °K for , 

2. 08 °K for CL and 2. 45 °K for NO while 0 is 3392 °K for N- , 

2 v 2 

2273 °K for 0 £ and 2738 °K for NO and 0 elec is 71584 °K for , 

11392 °K for and 174 °K for NO (ref. 27). In other words, in this 

range of temperatures, if we have or a mixture of the two, 

excluding nitric oxide, at equilibrium at some gas temperature T^ = 

S 

T^ ot = T v ^ = ^trais = etc * 9 S as > to S ooc * approximation, will have 

the same energy loss factor as a gas with the same vibrational temperature 
Tg = T^k k ut different rotational and translational temperatures, assum- 


ing that T . » 0 
& rot rot 


and that the vibrational and electronic excitation 
are at equilibrium at the same temperature (or that the contribution to the 
energy loss factor by electronic excitation is negligible). Since the 
characteristic temperature for electronic excitation for is high, the 

contribution to the electron energy balance due to electronic excitation is 
negligible in the temperature range of interest and, in this range, 


partners. The assumptions made in this derivation are: (1) ideal 
harmonic oscillator behavior of the diatomic molecules, (2) the Landau- 
Teller approximation for the collisional transition probabilities 
[Landau, L. ; and Teller, E. : Physik. Z. Sowjetunion 10, (1936) p. 34], 
(3) the existence of a Boltzmann distribution for both the population 
of the vibrational states and for the energy of the free electrons. 
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knowledge of the variation of 6 M with vibrational temperature is suffi- 

™ 2 

cient to determine the influence of gas temperature on 6^ (always 
provided ^ TQt is much greater than ® TOt )• hi other worL, vibrationally 
hot Nt is a 1 < aimply hot as far as 6^ is concerned. The same state- 
ments can also be made for oxygen, at least in the lower half of the 
temperature range of interest (and in the upper half dissociation begins to 
make the question academic). 

Therefore, for nitrogen, oxygen or any mixtures (that 

somehow are prevented from containing NO ) at equilibrium at any tempera- 
ture T in the range from approximately 300 °K to 10, 000 °K, the 
electron energy loss factor is the same as the electron energy loss factor 
of the same mixtures at a vibrational temperature T^ = T^ in the same 
range but not at equilibrium, provided the rotational temperature is at 
least 200 - 300 °K. We see from table II that this condition is met in our 
experiments and therefore T ^ , which turns out to be essentially equal 

to T ^ ~ in this work, is sufficient to characterize the effective 

e, E = 0 

temperature T^ of the nitrogen gas as far as the electron energy-loss 
factor is concerned. 


The situation is not necessarily the same for other gases. For 
example, with monatomic gases it follows from the definition of 6^ that 
the translational temperature should be sufficient to characterize the 
energy state or effective temperature T of the gas as far as 6^^ is 
cor cerned. 


EXPERIMENTAL RESULTS 

Experiments were carried out to measure the electron temperature 
as a function of the power input to the cell, in order to obtain dP/dT e in 
argon, nitrogen and nitrogen-oxygen mixtures preheated by the arcjet. 

In all these experiments the reflected rf power was negligible and there 
was no observable interaction between the Langmuir probe and the micro- 
wave cavity or the rf power (e. g. , the reflected power did not change when 
the probe was withdrawn from the cavity). The ohmic heating power 
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loss SI was also negligible compared to the illumination power per unit 
volume of plasma. 

As far as could be determined, the species concentrations at the 
probe were not at equilibrium with , or with T g g - q * or an 

easily identifiable excitation temperature of the gas T* , or with the 
electron temperature after the application of rf power, or with the 
stagnation temperature of the gas in the plenum chamber, T q , as 
determined from an average enthalpy derived from an arcjet energy 
balance. All evidence pointed to the existence of an arc temperature 
T arc > T q which was responsible for the species present originally in 
the arc chamber. 

Typically, in these experiments the particle concentrations were 
n a 10*^ electrons/m^ or slightly larger and n a 10^ atoms or 
molecules per in' 3 . The electron concentrations were obtained by the 
Langmuir probes and the neutral particle concentrations from the pressure 
in the tank and the static temperature in the jet. Note again, however, 
that in accord with equation (20) these measurements are not required to 
find the energy loss factor. They serve to define the range of experimental 
conditions and to verify that the constraints listed earlier are in fact 
satisfied. 


Experiments with Argon 


For the argon experiments, the appropriate gas temperature T^ 
c T for equations (4) or (8) is the temperature that describes the 

Kj S 

average energy (excluding electronic excitation with the associated tem- 
perature, ) of the atoms and is identical to the static temperature 

T of the gas. In a simple gas, consisting of one heavy component, 6^ 
is of course equal to the energy loss factor 6^ of that gas. Choosing the 
reference temperature T to be the gas temperature T , which is 

^ O 

identical with T in this case, we can use equations (7) or (9) to define 


6 k = 6 eff = 6 A for ar g° n - 


Then we can use the appropriate expression 


among equations (17) to (20) to determine 5^ from the experimental data. 
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When this experiment was carried out with argon under conditions 
such (ref. 28) that the electron temperature at the probe with zero applied 
electric field, T r? _ n > was equal to the excitation temperature T" 1 = 

C , JcL» — U 

T e ie C the argon atoms and was in the range 2100 K < T g ^ ^ < 

o ^ 

3400 K, it was found, by using tie second half of equation (20) in a 

simplified form given by 


6 eff " 2.86 x 10‘ 6 (^) (£?-) = 0- 312 Tjr-^jr— (43) 

r e e e, E = 0 

that the computed energy-loss factor was approximately 13% higher, at 
_ 5 

6 e ^ = 3. 10 X 10 , than the elastic energy-loss factor for argon, 

6 a . = 2. 73 X 10’^. Since the enthalpy of the argon in these experiments 

^ 6 

was in the range of 10 joules/kg to 2 X 10 joules/ kg, while the static 

temperature T = T = T of the gas was in the range of 100 to 200 °K, 

S * 

it appears unlikely that this discrepancy was due to radiation or other 
losses. It is much more likely that the electric field in the waveguide 
(or the illumination power applied to the plasma) is overestimated by this 
amount. The results of the measurement of the waveguide electric field, 
discussed in the section on experimental apparatus, appear to confirm 
this view. 

Typical results with argon at these low enthalpies are summarized 
in table III. The stagnation temperature of the argon was in the range of 
2000 to 3500 °K. The results of table III are obtained by making use of 
equation (43). There is clear evidence (ref. 1) that at higher electron 
temperatures T v - n » the energy-loss factor for argon as c ;• nputed 
from equation (43) is larger than the elastic energy-loss factor. This 
result indicates that some inelastic process or processes serve as a sink 
of the electron energy in argon and that these processes significantly 
influence the determination of 6 ££ by equation (43) due to the smallness 
of 6 ^ and of the term (3/2) kn g v t ^ e ff equation (34) for monatomic 
gases. 

To illustrate the significance of radiation energy losses note that 
22 3 

at 2. 7 X 10 argo i atoms per m (corresponding to 370 microns tank 
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pressure and a free jet static temperature of 100 K) and an electron 
temperature of 10, 500 °K, the total line radiation intensity from an argon 
plasma is approximately 3. 0 watts/cm (ref. 29) if the population of the 
radiating states is in equilibrium with the electron temperature. At an 
illumination power of 1 watt, this radiation loss is about 1000 times 
larger than the ohmic heating given by equation (41) for argon. However, 
to obtain a truer measure of the significance of radiation energy losses, 
we must compare the term (3/2) kn^ ^ e ff * with units of W/°K - m^ , 
to the term 3ft/9T e in equation (34). To make this comparison we assume 
that the line radiation intensity from argon changes with temperature at 
this density at the same rate as the line radiation from oxygen or nitrogen 


_ T* 


Then use can be made of the data of Allen 


atoms and that T = T* 

ti o 22 

(ref. 16) to compare these two terms. At T = 4900 K, n = 2, 7 X 10 

3 - 20 2 ^ a 8-1 
argon atoms per rn and *2X10” m , we obtain = 2 X 10 sec 

Then for n = 7 X 10 *^ an d 5 = 2. 7 X 10”^ we obtain (3/2) kn v.6 ,, = 

<j C a A q a Cl CXX 

7.8X10 W/°K - m or 7.8X10^ W/°K - cm . From the data of Allen 

we obtain at T V ~ 4900 °K and at the same number density of argon atoms, 

AR/ AT* = R 490q - R 3900 /1000 * R^^/AOOO « 40 X 10 ' 9 W/°K - cm' 

(ref. 16). Therefore, the term 3R/8T e appears to be much larger than 

(3/2) kn g v t 6 g ££ at 4900 °K and cannot be neglected. Naturally, the 

excitation temperatures T* (which are usually equal to T g ^ _ q under 

the conditions of the experiment) are not quite as nigh in the results of 

table III. However, the departure of the argon gas in the jet from local 

thermodynamic equilibrium is so strong, as evidenced by the magnitude 

of T_ it _ a > that only actual measurements of the radiation losses from 
e , hj - U m 

the gas can settle the question of the influence of 3R/9T e and instill 
confidence in the precision of the results of the determination of 6 ^ . 


Careful estimates of the rate of ionization and recombination 
within the waveguide under the conditions of the experiment seem to ex- 
clude the possibility of significant errors in the measurements due to these 
reactions. However, at this time we can only speculate on the influence 
of the concentration and rates of excitation and de -excitation of metastable 
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argon atoms on the magnitude of the term 
measurements. 


2 s.e. and its effect on the 
< x 1 


An alternative approach would be to include the effect of radiation 
losses and reaction processes in the quantity 6 ££ determined by these 
experiments as a function of and by summing the terms 

(3/2) k n e 6 e f £ and 3 3?/3T e and expressing them in terms of an equiva- 
lent 3 e ££ , 8 , defined by 




( 


2 

3 k n 

e 



3 $ 
cFF 


(44) 


This is a highly artificial approach since it masks many important effects. 
For example one cannot assess from it the influence of geometry and the 
relative importance of radiation losses and reaction processes. Therefore 
results obtained in this way are of doubtful predictive value. 


Experiments with Nitrogen 


When the experiment was carried out with nitrogen under such 
conditions that, at the probe location, the nitrogen was molecular and the 
electron temperature with zero applied electric field, T g - q ’ was 
always equal to an excitation temperature T r that could be closely identi- 


fied to the vibrational temperature, T 


vib 


of the nitrogen molecules 


(see table II), it was found that the temperature T g e - 0 = ^ = ^vib 

was always much higher than the static temperature computed by an 
isentropic expansion from the conditions of the plenum chamber (stagna- 
tion pressure P q « 0. 5 to 1.0 atm and stagnation temperature T q = 1500 
3500 °K) to the vacuum tank pressure (200 to 1000 microns). In addition, 
when the vibrational relaxation time was longer than the time of flight of 
the molecules from the plenum chamber to the probe, the vibrational 
excitation temperature T^, for the nitrogen molecules at the probe was 
approximately equal to or higher than the stagnation temperature T q in 
the plenum chamber (but always less than the excitation temperature in 
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the arc). Thus in the absence of an electric field, the electron temperature 

is in equilibrium with the vibrational excitation temperature of the 

molecular nitrogen T and we can write T ~ n = T £ T # There - 
e vib e, E = 0 vib o 

fore, we can write T = T = T _ n = T ., in equations (5) and (17) 

n g e , E = 0 vib 

and the appropriate temperature that characterizes the energy state of the 
gas is T v ^ = T e £_q a s far as the energy-loss factor is concerned. 

The coupling and equilibration of the free -electron and vibrational 

temperatures under essentially similar conditions (but in the presence of 
an excess of argon) has also been established by Hurle and Russo (ref. 30). 

_ 5 

Note that 5^ ^ = 3. 9 X 10 while Crompton and Sutton (ref. 18) 

give, for 6„ T , at T* = 288 °K, 6„ = 3. 33 X 10" 4 at T = 893 °K, 

N 2 ' g N 2 a e 

=2.90X10 at T =3140°K, 6„ T = 3. 28 X 10" at T = 5260 °K 

n 2 4 e 9 N 2 e 

and 6^ = 8. 20 X 10 at T^ = 9060 °K. As a consequence of the fact that 
at Tg = 288 °K the energy-loss factor for N ? is nearly constant over 
the temperature range 893 °K < T^ < 5260 °K we can use equation (20) 
in a simplified form given by (see equation (43)) 


6 ,, = 0. 287 

eff 


e, E = 0 


where the constant 0. 287, corresponding to Z = 387 ft , is the result of a 
compromise that partly compensates for the high values of 6^ determined 
when an impedance of 420 ft is used. The lower value of Z is also used 
in all determinations of 6^ for mixtures reported here. 

Typical results of the variation of the measured electron tempera- 
ture with rf power in nitrogen at 0. 9 mm Hg are shown in table IV. For 
these runs T q = 3230 °K, as obtained from an average enthalpy derived 
from an arcjet energy balance. Typical values of 6^ as a function of 
gas temperature are shown in table V. The results of five or more runs 
are averaged for each T^ , with typical spreads as in table IV. Equation 
(45) was used to obtain all data in tables IV and V, although the last four 
entries of table V are slightly beyond the range where equation (45) was 
intended to be used. 
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As far as could be determined the species concentrations at the 
probe were not at equilibrium with ^ - g " ^vib at a m uch lower 
temperature (see also reference 12). Therefore, the gas at the probe was 
almost pure N 2 and the values of 6 ^ tabulated in tables IV and V 


should correspond to 6 


N 2 


However, although "resonances 11 in 6^ with 


T cannot be excluded, it is felt that the variations in the tabulated results 
§ 

are due to impurities in the gas (welder's grade bottled gases, water 

pumped, with large '•■ariations in impurities) and/or leaks in the apparatus. 

It was observed, for example, that the values of " ^ obtained in these 

experiments changed (a) quite frequently, by small amounts, when the 

exhausted gas cylinders were replaced by new ones and (b) always v/hen 

oxygen was deliberately leaked into the gas feed line in small quantities 

(approximately 0. 5% ). Air leaks were found in the apparatus sometime 

after the two series of runs corresponding to = 1750 and 3480 °K of 

table V were carried out. The observed increase of 6 rr = 6^ at 

eff N 2 

T = 5950 K may be due to increased elastic losses due to the rather high 

^ o 

electron temperature. The observed decrease of 6^^ at T^ = 6100 K 

may be due to an increase of the atom concentration. Whether indeed an 
increase of 150 °K in T^ can change the atom concentration sufficiently 
to cause the observed shift in the value of 5^ is still a matter of con- 
jecture. On the other hand, the temperatures in these last few entries of 
table V are beyond the range where equation (45) was intended to be used 
and therefore this behavior is not surprising. 


Typically in these experiments with nitrogen we would expect the 
quantity (3/2) kn^v^^ to be larger by one or two orders of magnitude 

than the value it has in the argon experiments. Thus at T = 3 600 °K, for 

22 3 16 3 e 

n = 10 N y molecules m , n =7X10 electrons/m and Q 6 ]st = 

9X10’ m^, we obtain (3/2) ^ n e v ^^ e ff = 2. 4 X 10 1 watts/ K - m . 

Since the main contribution to the energy radiated from a nitrogen plasma 

at these temperatures and pressures comes from (1 +) and atomic line 

radiation (i. e. , from electronically excited states) and since the radiative 

lifetimes for spontaneous emission for these transitions is short (of the 

order of 10 y sec or less) and the population of these excited states is at 
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equilibrium at a temperature much l~wer than T — , the radiation 

, e * = U 

term 8R/3T e and the term 8 (S s j e ^)/9T e can be neglected for all 
illumination power increments more than 0. 2 watts. 


As a consequence of the discussion elsewhere in this report, 

energy-loss factor for hot nitrogen at equilibrium at a temperature 

is expected to be the same as the energy-loss factor at T = T ^ 

sy g e, E = 

tabulated in table V. 


the 

T 

g 

0 


Experiments with Nitrogen-Oxygen Mixtures 


Typical values of 6 ££ for N^-O^ mixtures as a function of the 
mole fraction of oxygen, = [ O^]/ ( [ O^] + [N^] ) , where the quantities 

in brackets denote particle concentrations in number/m\ are shown in 
table VI and figure 3. These results were obtained by setting the arc jet 
power controls so as to attain a stagnation enthalpy of approximately 

£ L 

3 X 10 to 4 X 10 joules/kg with pure N ? (corresponding to T = 2500 - 
3200 °K) and then introducing an increasing amount of in the secondary 
flow to the plenum chamber while the total flow rate was kept approxi- 
mately constant by simultaneously decreasing the flow rate of secondary 
. The primary flow (i. e. , the flow through the arc chamber) was pure 
and was kept constant. Mixing of the two streams in the plenum 
chamber was assured by a long-enough residence time (of the order of 
1 m sec) and the evidence obtained by color photography, radial surveys of 
the ion and electron currents and the electron temperature in the plasmajet 
seemed to confirm that mixing of the two gas streams (i. e. , the hot 
primary ^ and the cold secondary mixture of and O ^ ) was perfect 

(see Appendix C). 


What is immediately apparent from the results of table VI and 
figure 3 is that even a small amount (~4%) of helps to de-excite the 
vibrationally hot pure ^ from an average = T g p _ q = 

345 0 °K (for the four pure nitrogen runs) to an average of approxi- 

mately 1950 °K (for the four mixture runs at * 0, 04). It appears that 
the degree of de- excitation increases and the effective gas temperature 
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decreases as the oxygen content increases but not by as large an 
amount. For example, from the four runs at an average » 0. 32 , we 
obtain T ~ 1550 °K. Therefore, the mixture data presented in table VI 

S Q 

and figure 3 correspond to an average gas temperature T of 1700 K. 

^ o 

For this reason, the pure N., value of 6.. at T = 1700 K should be 

2 N2 g 

taken as the value corresponding to the mixture data of table VI and 

figure 3. Reference 1 gives = 2. 9 X 10"^ at 1700 °K, We note that 

these values of appear reasonable. For example, for air (M^ ~ 0. 21) 

at room temperature, Ginzburg and Gurevich (ref. 7) obtain 6 . = 

« -2 air 

1. 2 X 10 at T = 2000 °K and 6 . = 1. 6 X 10 at T = 3000 °K, while 

e air ^ e 

Hake and Phelps (ref. 31) report 1. 7 X 10~ for the energy-loss factor at a 

"characteristic energy" of 0, 2 ev. These values compare well with 
_ 3 

6 a - r = 3X10 obtained from figure 3 for air at a temperature in the 
vicinity of 1700 °K. 


An accurate estimate of the term 8$/8T in N->-CL mixtures 

e 2 2 

cannot be made since the radiated power is a strong function of the con- 
centration of nitric oxide present in the mixture and there are no reliable 
estimates of [ NO] . On the other hand, the magnitude of the term 
(3/2) k n e v ^ e ff is even larger in these mixtures than it is in pure 
nitrogen and it is possible that it remains larger than 8$/8T e for all 
mixture ratios. 


DISCUSSION OF RESULTS AND CONCLUSIONS 

There is not much else that needs to be said concerning the results 
for argon and nitrogen presented in tables II and V. If anything, these 
results are surprisingly close to the expected values for argon and appear 
reasonable for nitrogen in the sense that they are within - 10% of the re- 
ported room temperature values at low gas temperature (T < 2000 °K) 

O 

where appreciable vibrational excitation has not yet taken place and they 
show an increase of for hotter nitrogen as compared to nitrogen at 

288 °K, a trend expected from theoretical considerations [ see, for 
example, equation (42)]. 
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Although the reproducibility and accuracy of the results for N^-O^ 
mixtures plotted in figure 3 could be considered satisfactory, there are 
two disturbing factors: When we use the consistent terminology of 
Demetriades and Argyropoulos (ref, 3) for the definition of collision fre- 
quencies t and cross-sections Q for collisions between electrons 
^ e, k eK 

and heavy particles, the mixture rule of reference i (assuming only N~ 

- 4 ^ 

and C >2 present), and the values of 6^ = 3X10 , Rq = = 

2. 5 and 2. 0, we derive for the energy-loss factor 6 of pure oxygen 
in the vicinity of 1700 °K the values listed in the last two columns of 
table VI. These values of 6^^ (2 X 10’^ to 3 X 10"^) are a little high. 

Thus if 5 q = 6. 7 X 10 ^ when T^ ~ 2000 °K and the gas is at room tem- 
perature (ref, 7), and if the pattern set in the nitrogen experiments (ref. 1) 
is followed by , one would not have expected 5 q^ to be much higher 
at 1700 °K than at room temperature. A clue as to the possible source of 
this anomaly is provided in table VII, where we see that the values derived 
for Sq^ by the mixture rule are dependent on the oxygen content. The 
source of this dependence appears to be the effect on the computation of 

of the presence of nitric oxide, NO, and atomic oxygen, O, in the 
mixture. In other words the assumption that or are the only 

species present breaks down. For example, we can show that at the lower 
oxygen concentrations, the assumption of a small concentration of NO 
is sufficient to decrease the computed value of for pure O^ by a 

factor of 2 or 3, because of the largeness of 6 ^q (ref. 32), even though 
the contribution of NO to the observed 6^ for the mixture is small. 

At the higher oxygen content, the assumption of a small concentration of 
atomic oxygen is sufficient to increase the computed value of 6 q by a 
smaller but still significant amount. It is quite likely that both species are 
present in small but significant amounts as far as the computation of 
6 q^ by the mixture rule is concerned. Electron attachment to O^ may 
also make a significant contribution to the observed values of 6 ££ . These 
are typical of the difficulties caused by the proliferation of variables due 
to chemical changes. Nevertheless, good accuracy can be expected in the 

observed 5 data for mixtures since the concentrations of both NO 
— - elt 

and O are probably small. 
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Careful determinations of the composition in the plasmajet or 
operation with other gases (such as argon-nitrogen and argon-oxygen 
mixtures instead of N^-O ^ ) are expected to remove this source of 
uncertainty in the data. 

The other disturbing factor is the rather unusual spread in the 
data at higher oxygen concentrations, especially when the intrinsic sim- 
plicity and accuracy of this method for the determination of is taken 

into consideration. Several sources appear to be contributing to this 
spread. Some of these are: (a) oxidation and surface contamination or 
aging of the tungsten probes; (b) noise in the arejet; (c) the increase in the 
number of variables especially in such mixtures as N^-O^ that makes it 
difficult to maintain one quantity constant while varying others*; and 
(d) the human factor in data reduction (especially of the Langmuir probe data 
where a certain minimum level of practice is required). All these sources 
of uncertainty, including the influence of 9<&/9T^ and especially the effect 
of radiation losses, need to be systematically attacked since it is felt that 
even small variations in the value of 6 ^ with gas temperature can have 
a very significant effect on the performance and stability of J X B and 
other devices (refs. 2 f 33). 

On the whole, the mixture energy-loss factors determined directly 

by this initial experiment and presented here are considered accurate within 
± 1/2 

a factor of (2) ' or better at low concentrations. At higher 0£ 

concentrations (M f > 0. 3) the additional complication of possibly significant 
departures from a Maxwellian electron energy distribution may have con- 
tributed to the spread of the data. Under such conditions the lower rf 
power runs should be more reliable (since they perturb the initial 
Maxwellian least). 

* For example, in these mixtures 6^ is a function not only of and 

T but also of [0 9 ]/ [N 9 ] , the stagnation enthalpy and pressure in the 
g c L 

arc chamber and the pressure in the vacuum tank - as a result even the 
job of presenting the data becomes very complicated. 
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The values for the energy-loss factor for heated pure oxygen derived 
by processing the mixture data are not as reliable but it would be sur- 
prising if the exact values do not fall within the spread of data, - ( 2 )* * , 
presented here. 

It is believed that the main sources of error associated with the 
equipment are: (a) drift of the Langmuir probe characteristics due to long- 
term fluctuations in the arcjet current- voltage characteristics caused by 
the wandering of the arc attachment spots on the anode and cathode and 
(b) noise due to intermittent arc operation or rotation as well as 60-Hz 
ground loops. The extent to which these sources of error affect the 
reproducibility and accuracy of the measurements is still a matter of 
conjecture. 

The fact remains that the experimental results follow very closely 
the expected values of the energy -loss factor for three different gases over 
two and one -half orders of magnitude. Moreover, this method for deter- 
mining the electron energy-loss factor, and therefore also the collision 
crons- sections for energy transfer between slow electrons and heavy 
particles, possesses two very significant and favorable features: (1) all 
important measurements are dc and (2) since the energy-loss factor has 
been uncoupled from (a) the collision frequency, or (b) the drift velocity 
and the diffusion coefficient, it can be determined independently of collision 
cross-sections for momentum transfer and particle densities. 
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APPENDIX A 


DIMENSIONS OF WAVEGUIDE AND ESTIMATE OF WAVE IMPEDANCE 


In choosing the dimensions of the waveguide several factors should 
be taken into account. These are: 

(1) The operating dominant mode of electromagnetic wave propaga- 
tion in the waveguide, TE^q , will excite higher modes because of the 
presence of obstacles (quartz tube containing the plasma, probe, etc. ) 
inside the waveguide. In the steady state not all of the higher modes carry 
real power because, depending on the dimensions of the guide and the 
operating frequency, some of them can be made evanescent. The TE^q 
scattered mode will enter in the series expansion for the scattered field 
(in terms of the waveguide modes) with the largest coefficient and can be 
termed the reflected wave as though only the TE^q existed. Then an 
estimate of the coefficient of the reflect 1 wave can give information about 
the coefficients of the higher modes that can be excited. 

(2) The holes on the sidewalls through which the plasmajet flows 
should not appreciably obstruct the flow of currents which sustain the pro- 
pagating mode. 

(3) The holes in the sidewalls should be large enough so that they 

* For a more detailed review of the subject of electromagnetic wave pro- 
pagation, waveguides and related measurements see (1) Ginzton, Edward 
L. : Microwave Measurements, McGraw-Hill Book Company, Inc. , New 
York (1957); (2) King, Donald D. : Measurements at Centimeter Wave- 
length, D. Van Nostrand Company, Inc. , New York (1952); (3) Reich, 
Herbert J. ; Ordung, Philip F. ; Krauss, Herbert L. : and Skalnik, John G. : 
Microwave Theory and Techniques, D. Van Nostrand Company, Inc. , New 
York (1953); (4) Wind, Moe; and Rapaport, Harold (eds. ): Handbook of 
Microwave Measurements, Microwave Research Institute, Polytechnic 
Institute of Brooklyn (1955); (5) Sommerfeld, Arnold: Electrodynamics, 
Academic Press, Inc. , New York (1964); (6) Goubau, George: Electro- 
magnetic Waveguides and Cavities, Pergamon Press, New York (1961). 
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do not appreciably obstruct the flow of the free jet of plasma. 

(4) The characteristic time of flow of the pla3ma across the wave- 
guide should be much greater than the relaxation time of the electron tem- 
perature to ensure an almost instantaneous steady- state response of the 
electron temperature to the local field so that the formula giving T - T 
in terms of the electric field be applicable. 

Let us consider the influence of each one of these design con- 
straints in selecting the dimensions of the waveguide. 

The magnitude of the wave vector of the TE nm mode in a rectan- 
gular waveguide is given by 


k 



) 


2 -. 1/2 


(A- 1 ) 


where a , b are the width and height of the waveguide respectively (see 
Fig. 1, Appendix A), n , m are integers not both zero, = w/c = 

51. 34 m * and o - 1. 54 X 10^ rad sec at f = 2. 45 GHz. 

A propagating mode should satisfy the inequality 

k^ > (mr/a)^ + (rmr/b)^ (A- 2) 

therefore, we may choose a , b such that this inequality is not satisfied 
for a given n , m and consequently any given TE^ mode will be 
evanescent, carrying no real power. For example, if we want only the 
TEjq mode to propagate down the waveguide we should choose 2b < \ Q = 

0. 122 m < 2a . The condition 2 a > \ ensures that TE^q propagates but 
does not exclude other modes from propagating. However, because of 
factors (2), (3) and (4) we cannot choose a , b such that only the dominant 
mode TEjq should propagate. Thus factors (2) and (3) require that 
b S: 4" = 0. 1016 m and factor (4) that a S 6" = 0. 1524 m . If we, there- 
fore, choose a = 6" = 0. 1524 m, and b = 4" = 0. 1016 m, we can easily 
show that the following modes can propagate and consequently carry real 
power: 
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TE 


20 ; 


TE 


01 


TE 


11 


(We do not have to worry about TM modes because we know that even if 
they get excited by the TE^ mode they will have a negligible coefficient 
simply because the exciting mode is a TE^q mode and not a TMjj mode. ) 

We now make an estimate of the coefficient of the reflected wave 
( = scattered TE^ ). 


The tube containing the plasma is made of quartz with an index of 
refraction approximately 2 in the microwave frequency range of the experi- 
ment. The index of refraction of the plasma in a good approximation is 
given by r\ = [ 1 - (w^/co 4 *)] l/ ^ 2 . 

2 2 

For a maximum gj equal to 0. 25 go we obtain r\ ~ 0. 8 67. The 
P 

plasma column will have a diameter of about 2. 5" or 0, 0625 m therefore 
the ratio of the power reflected to the incident power for a plasma slab of 
thickness d - 2" or 0. 05 m is 


where 


2s [ 1 - cos (2 k d Ti)] 

p = U 1 

1 + s - 2s cos (2 k d r\) 


s 



k 



v'ir/ a) 2 ] 


1/2 


46 m~ * 


(A- 3) 


Therefore the maximum power reflected from the plasma p^ is 1. 7 per- 
cent of the incident power. For a quartz thickness d ~ 0. 001 m the 
reflectivity of the quartz (slab) pQ will be ~ 0. 58 percent. The resulting 
p is of the order of the sum p + Pq because the index of refraction of 
the plasma is close to unity. 


The above rough estimate of the coefficient of the reflected field 
and power assures that the coefficient of the higher modes will be small 
compared to the TE^q mode so that the power carried by the non- 
evanescsnt scattered modes will be negligible compared to the power car- 
ried by the incident dominant mode TE^q . This last conclusion is further 
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supported by the fact that the wave impedance of a mode increases with the 
order of the mode (see below). 


Wave Impedance 


When 86 « 0, as we have already explained in the text, follow- 
ing equation (19), the electron energy-loss factor is given by 

-1 


eff 


T7 2 \7T 

3 km co r 

e 


n ■ (£) 


(A-4) 


Where Z = impedance of the gap in the waveguide, P = power and A = 
the cross-sectional area of the waveguide. Now, in general, a propagating 
TE nm mode, when n , m + 0 , carries a power given by 


nm 


ab 

T 


{ |E | 2 + | E | 

1 OX 1 I nv I 


"TZ 


oy 


h, nm 


or, when m = 0 and n + 0 , 


(A- 5) 


with 


where 


nO 


ab 
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oy 
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kZ 


J h, nm 
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ii 
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TyTJT 
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i 
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= = k 


(A- 6) 


(A-7) 


where p , e , r\ , are the permeability, permitivity and index of 
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i /2 

refraction [ r\ = (jjl e /jx e ) 7 ] respectively of the medium by which the 

guide is filled. 


For vacuum and the dominant mode we obtain Z, 4 = 1. 09 Z 

1/2 ^ ^ ^ 

where Z = (u /e ) 7 = 377 £2 . For the plasma in question with max gj 

o *0 o r p 


(worst case) we obtain Z, .- = 1.3 Z 

h, 10 o 

should be such that 


Therefore the actual Z 


h, 10 


l ' WZ o < Z h, 10 < 1 * 3 Z o (A - 8 > 

and of course it will be closer to 1. 09 Z than to 1. 3 Z . (We do not 

o o 

have to consider the effect of the quartz tube since 2-n-k = \ » d where 

d is the quartz thickness). 

Now from equation (A- 6) the amplitude of the electric field will be 
given by 



o 


4 P 

IF * 


Z h, 10 


(A-9) 


and Z in equation (A-4) is properly given by Z^ ^ while the area in 

equation (A-4) is given by A = ab/4 . The average electric field E^ will 
2 2 2 

be given by 2E = 2 (E ) = E q . The origin of the factor 4 in equation 
(A-9) can be interpreted to be the result of two averagings over time and 
space . As P in equation (A-4) we may use the power transmitted through 
the guide. This power will be measured by the apparatus we have designed 
and will be given to good approximation, at the center of the waveguide 
where T^ is measured, as the power input to the waveguide plus the power 
output divided by two (this accounts for connector and waveguide losses). 
Thus P usually amounts to 80% of the input or nominal power. 

The power losses in the plasma have been shown to be negligible 
under the conditions of the experiment. 

2 2 

Numerical values of the grouping [e /(Ikm^u )] (4Z/A^) that 
appears in equations (19), (20) and (A-4) as a function of the impedance Z at 
f = 2. 45 GHz are given in the following tabulation: 
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Z = Z = 377 n, 
o 

Z = 387 ft , 

Z = 400 £2 , 

Z = 1. 09 Z = 411 £2, 
o 

Z = 1 . 3 Z = 490 n. 


6 , f = 0. 280 (AP/AT ) 
eff e' 

6 eff = °- 287 (AP / AT e } 
6 eff = 0. 297 (AP/AT e ) 

6 eff = 0. 305 (AP/AT e ) 
6 eff = °* 363 ( AP / AT e > 


Numerical values of some other groupings of interest are given in 
the following tabulation for f = 2. 45 GHz (to = 1. 5393 X lO 1 ^ rad/sec) and 
mks units: 

— = 5. 947 X 10 -29 , = 0. 4829 X 10 23 

2m co 
e 


= 0. 2871 X 10 


-5 


3km to 
e 


3 k 6 


= 5. 645 X 10' 


28 


-2 2 

while for Z = 387 ohms and A = 0. 1016 m X 0. 1524 m = 1. 548 X 10 m 
we obtain 


E 2 = (4Z/A r )P = 1.0X10 5 P 

where P is in watts. 
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APPENDIX B 


SENSITIVITY OF THE MEASUREMENTS 


At 2. 45 GHz (w = 1.54X 10^ rads/ sec) the quantity 2e / (3 km to 2 ) 


is equal to 4 X 10 15 , (3 X 2. 4 X 10 20 ) = 5. 74 X 10' 6 , mks units. 


Therefore 


if T g = f (E q ) is measured we have from equation (20) in mks units 


6 = 


5. 74 X 10‘ 6 E 


,dT 


-1 




Let us assume that at a given gas temperature T a plot of T is obtained 

§ ® 

as a function of E q and that for our purposes E = 0. 707 E q . Then we can 

read the corresponding slope dT e /dE Q and E q at any T e and obtain 6 . 

We have said that the quantity 6 should not vary more than a few percent 

within narrow ranges of T . Therefore at a given T the quantity 

e - 1 8 

E • dE /dT or E (dT /dE )" should be constant at approximately 

5 o 

2X10 6 in mks units. For example, if, at T = 4000 K and a given 

T , the electric field is E =10 volts/cm =10^ volts/m and dE /dT = 

0.5 volts/(m-°K) or dT e /dE o = 2 °K-m/vclt , we obtain 6 = 2. 8 X 10* . 

Then also at T g = 3000 °K (and the same T^ ) and at E q = 5 volts/cm, 

we must read approximately dE Q /dT e ® 1 or dT g /dE o * 1 in mks units. 


If T can be 
e 


The slope of T g vs. E q will be proportional to E q . 

measured within ±100 °K, a range of at least 1000 °K must be available 
for T g variation to ensure good accuracy and the required change of E q 
over the interval must range, conservatively, between 0. 1 volt/m per °K 
and 1 volt/m per °K. Since it will be possible to vary E q by hundreds 
of volts/m (1 volt/cm = 100 volts/m) we should be able to obtain changes 
of one thousand or more °K in electron temperature and therefore it seems 
possible to obtain values of 6 with good precision. Operation at higher 
T g and E q will enhance this precision if the error in measuring T g 
remains the same. 


In this Appendix, we omit the subscript "eff" from 6, for simplicity. 
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From equation (17), with v. » oo , we obtain 


6 


2 

e 

3 km oo 
e 



T) 


Now let us consider the error, A6 , in 6 , resulting from the uncertainties 
AE q and AX (X = T g - T n ) when 6 is determined from a single measure- 
ment of E q and X . The error A6 will be given by the familiar ex- 
pression 


now 


hence 


A6 = 


96 
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96 
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96 
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This gives the relative error in 6 which is 


A6 _ AX , o 

~T~ ~ IT E 

o 

Since AX = | AT | + | AT |, if we can assume that AT = 100 °K, 

6 a n © 

AT = 50 °K and X s T - T = 10 K while AE =100 volts/m and E = 

^ e n o o 

10 volts/m, we obtain, for the relative error in 6 


100 X 


A6 

TT 


100 X 


( 100 + 50 

l 10 0 0" ■ 


2 X 100 \ 

“ioo'cr ) 


35% 
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It is important to realize that this error is the maximum error that 
can occur at these values of E q and X and when only these two measure- 
ments are made. Even for such single measurements there is a probability 

that the errors in T , T and E will partially cancel. The situation is 

e n o r 7 

even more favorable for the experiments described in the text because many 
measurements of differing E q and have been made in determining the 

value of 6 corresponding to a given T n (i. e. , in determining the slope 
dT e /dE Q ). When 6 is obtained from such multiple measurements then the 
error is greatly reduced by averaging effects. 

An evaluation of the results presented in the literature on measure- 
ments of 6 = f(T g ) appears to indicate that other methods for evaluating 
6 (by swarm techniques) are less accurate than the method described here. 
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APPENDIX C 


EXPERIMENTAL PROBE CHARACTERISTICS; SURVEYS OF ELECTRON 
AND ION CURRENTS AND ELECTRON TEMPERATURES; TYPICAL 
ROTATIONAL AND VIBRATIONAL SPECTRA 

Figures 1C and 2C of this Appendix are traces of the probe current- 
voltage characteristics obtained for nitrogen on the x-y recorder with and 
without rf illumination power. In figure 1C the characteristics were traced 
twice and in figure 2C five separate times for each recorder sensitivity 
setting. This procedure of tracing the probe characteristics more than once 
on each recording sheet was established and carried out routinely in all 
runs in order to provide a visual indication of drift or excessive noise. 

Figures 3C and 4C of this Appendix represent a radial and axial 
survey, respectively, of the electron and ion currents in the plasmajet with 
the Pt-Ir gas-cooled probe. The radial survey, figure 3C, was carried 
out by traversing the probe on the z-axis, i. e. , perpendicular to the 
plasmajet axis (x-axis), at a distance of approximately 0. 35 m downstream 
from the arcjet nozzle at a fixed operating condition (i. e. , constant mass 
flow rate and power) corresponding to a stagnation temperature of approxi- 
mately 2500 °K. The uncertainty in locating the jet axis is less than 
±0. 05". Note that the probe was traversed from left to right and then 
returned from right to left (i, e. , it started 2" to the left of the plasmajet 
axis as we look down the axis of the tank from outside the tank door and 
continued traversing past the jet axis or centerline to 2" to the right of the 
plasmajet axis and then returned back to its original position 2" to the left 
of the plasmajet). Note the perfect azimuthal symmetry of the electron and 
ion current distributions. Since this symmetry of the jet was obtained 
when 50% of the flow rate was injected into the plenum chamber downstream 
of the arc, it is concluded that the primary and secondary gases mix 
perfectly in the plenum chamber. 

In the axial survey, figure 4C, the sharp peaks in electron current 
correspond to passage of the probe through the apex of the shock diamonds 
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in the supersonic plasmajet. Note that the peaks in ion current lag the 
peaks in electron current. This axial survey was obtained by placing the 
probe at the exit plane of the nozzle and withdrawing it along the jet axis. 
This procedure was followed for both the electron and ion current traces. 

The probe location was checked repeatedly during the production of both 
traces and there was perfect spatial coincidence, i. e. , there can be no 
question of the ion-current trace being displaced with respect to the electron 
current trace. The coincidence of the bottom of the valleys in both traces 
confirms the spatial coincidence of the two traces. No drifts were observed 
in the arcjet characteristics during the recording of these traces. A total 
of eight clearly-detectable probe-current peaks were obtained when the 
probe was traversed up to a distance of 35" from the nozzle exit plane, 
corresponding to an equal number of shock diamonds. 

Clearly, the results of the axial survey indicate that the probe 
should not be located near or at the apex of a shock diamond since large 
concentration and/or electron temperature gradients occur there and the 
conditions leading to equation (16) may be violated. 

Figures 5C and 6C are the results of a radial and axial electron 
temperature distribution survey, respectively. 

Comparing the results of figures 3C and 5C we see that although 
the electron temperature remains very nearly constant across the jet, 

(from + i. 0" to - 1, 00") the ion current goes from essentially zero to its 
maximum of 0. 375 mv in the same interval (2,5 X 0,81 * 2. 0"). At the 
same time, the electron current decreases less rapidly than the ion 
current. Thus the ion current reaches 10% of its peak centerline value 
within 0. 8" from the centerline or jet axis whereas the electron current 
reaches 10% of its peak value at approximately 1. 12" from the centerline. 
The implication of these observations is that a low-n g high-T g sheath 
surrounds the jet of plasma. 

Another interesting observation concerns the increased noise level 
of the electron current signal as the probe traverses the jet boundary 
and/or the sheath that becomes apparent in figure 3C. This explains why 


some early runs that were made with a probe that was not carefully 
centered had to be scrapped because of too much noise. 

The lack of a clear-cut fine structure or symmetry in the T g dis- 
tribution plotted in figure 5C was somewhat of a disappointment which was 
only partially compensated by the excellent linearity over 2. 5 to 3 decades 
of the probe data for the temperature points plotted in that figure and 
obtained with the high- sensitivity x-y recorder at 0. 1 mv/in. The plasma 
was pure nitrogen and the probe was the platinum-iridium gas-cooled probe. 
The uncertainty of each point plotted in figure 5C could not have been more 
than ± 5% and the spread in the data must be due to drift of the arcjet 
characteristics. 

Figure 6C depicts the probe characteristics obtained at different 
distances from the arcjet nozzle exit plane. The remarkable feature of 
this survey is the uniformity of the shape of the probe characteristic 
traces. The uniformity of the slope of these curves testifies to the uni- 
formity of the electron temperature with axial distance (x-axis), at least in 
the interval from 4 inches 0. 10 m) to 19 inches (~ 0. 48 m) from the 
nozzle exit. 

Note that if the electron temperature does indeed remain constant 
or nearly constant through the oblique shock diamonds of figure 4C the 
peaks in electron current must correspond to increased electron concen- 
tration and the peaks in ion current to increased ion concentration. How- 
ever, this might violate the charge neutrality condition. Therefore, the 
lag between the electron and ion current peaks must signify that some other 
significant relaxation effect is taking place - for example, the ion tempera - 
ture may be changing through the shock diamond apex. This phenomenon 
has not yet received a satisfactory explanation. 

The surveys of figures 2C through 6C were obtained in a free jet 
of nitrogen plasma corresponding to arc operating conditions similar to 
those for which T r> n K 3000 °K in table V of the text. The waveguide 

C , Jh = U 

and 2, 40-inch O. D. qu. rtz tube surrounding the jet were removed for 
these surveys. 
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Figure 7C depicts the probe characteristics obtained in an axial 
survey within the waveguide cavity (i. e. , with the quartz tube and wave- 
guide in place) for approximately the same arc jet operating conditions with 
rf power on and with rf power off. A clear-cut decrease of the slope of 
the curves obtained with rf power on as compared to the slope of the 
curves with rf power off and as the probe location moves toward the cavity 
center, is apparent. Surveys similar to figure 7C confirm the existence of 
an electron temperature distribution within the waveguide with rf power on 
of the form T g oc [E sin(2x/a)]^ . 

Figure 8C shows the effect of rf on the probe electron current. 

It was obtained by carrying out two traverses with the gas-cooled probe. 
Each traverse begins at the nozzle exit and continues to about 20" down- 
stream (i. e. , past the end of the waveguide). In one traverse the rf power 
is on and in the other it is off. The important results are that the two 
traces come together past the end of the waveguide (proving that the elec- 
tron temperature rise in the waveguide does not produce any persistent 
effects) and that the shock diamond peaks are enhanced with rf power. 

Figures 9C and 10C are typical results of the electron-beam 
excitation -spectroscopy diagnostic experiments obtained with the apparatus 
described in the text. Figure 9C shows the rotational band of (0, 0) at 
3914 A used to determine rotational temperature and figure 10C shows 
typical band profiles of the (0, 1) progression at 4278, 4236 'ad 4199 A, 
used to determine the vibrational temperature. 

Note that the interpretation of these intensity profiles and the probe 
characteristic traces leads to the conclusion that the vibrational and rota- 
tional energy states as well as the free electron temperature seem to have 

Maxwellian distributions so that T T . and T n are well- 

vib ’ rot e, E = 0 

defined. 
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TABLE I 


SUMMARY OF CONSTRAINTS 


Origin 


Expression 


Implication 


1, Maxwellian distribu- Equation (20) 

tion of electrons 
assumed 


2, Penetration of rf Equation (23) 

field into plasma 


3. Frequency of electric Equation (24) 

field much greater 
than collision frequency 


n /n £ 10 ’ 1 

e' k 

Sets an upper limit on 
n for a given n . Can 
be relaxed for small 
electric field. 

n < w 2 /324 it 2 
e 

Sets an upper limit on 
n^ for a given co . Can 
be relaxed if high elec- 
tron concentration core 
is only a few electron 
mean free paths thick. 

n 2 « i. 46 X 10“ 8 w 2 /Q 2 T 
a ea e 

Sets an upper limit on 
density of neutrals. 


4. Langmuir probe size 
smaller than electron 
mean free path 

5. Langmuir probe size 
greater than Debye 
length 


Equation (25) Sets an upper limit on 

density of neutrals. 

Equation (26) Sets a lower limit on 

electron density. 
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TABLE I 


(Concluded) 


Constraints 6-11 are associated with the experimental 
technique chosen to heat the electrons 


6. 

Adequate flow time 

Equation (28) 

Sets a lower limit on 


across waveguide 


density of neutrals 
(weaker than (2) and (3)). 

7. 

Adequate increase of 

Equation (30) 

Sets a lower limit on 


electron temperature 


amplitude of rf field. 

8. 

Assures validity of 

Equation (32) 

Sets an upper limit on 


equation (17) 


amplitude of rf field. 

9. 

Breakdown field not 


Same as (8) 


exceeded 



10. 

Binary collision 
theory to be valid 


Same as (8) 

11. 

Sufficient increase of 

Equation (36) 

Defines a permissible 


electron temperature 


range for increase of 


that radiation and re- 


electron temperature. 


action energy losses 
are small compared 
with other inelastic 
energy exchanges, 
but not so large that 
6 g ff changes much 
with increase in T 

e 
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TABLE II 


TYPICAL PLASMA TEMPERATURES IN °K FOR 
ARC HEATED PURE NITROGEN EXPANDED TO A 
PRESSURE OF APPROXIMATELY 400 MICRONS 


T (a) 
o 

T (b) 

e, E = 0 

T <c) 

1 vib 

T (d) 

rot 

irt _ t (®) 

tr 1 static 

1800 

1720 

1810 

275 

260 

1950 

1650 

1880 

295 

290 

2200 

2800 

2590 

305 

310 

2400 

2580 

2660 

320 

350 

2800 

3230 

2980 

385 

425 

3050 

3580 

3450 

420 

475 

3350 

3620 

3840 

540 

540 


(a) Stagnation temperature from average enthalpy from arcjet 
energy balance. 

(b) Electron temperature from Langmuir probe, electron beam off. 

(c) Electron beam spectroscopy. 

(d) Electron beam spectroscopy. 

(e) From isentropic expansion. 
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TABLE III 


TYPICAL ENERGY-LOSS FACTORS FOR ARGON AT A 
GAS TEMPERATURE BETWEEN 100 AND 200 °K 


Electron Temp. , 
T e, E = 0 ’ ° K 

rf Power, 
watts 

Electron Temp. , 
T g , °K (with rf ) 

Energy-loss 
Factor, 6 eff 

3395 

0. 217 

5660 

3. 00 X 10‘ 5 

3220 

0. 220 

5450 

3. 08 X 10‘ 5 

2420 

0. 285 

5330 

3. 06 X 10" 5 

2165 

0. 2675 

4655 

3. 35 X 10" 5 

2150 

0. 275 

4860 

3. 17 X 10" 5 

2515 

0. 2975 

5445 

3. 17 X 10‘ 5 

2335 

0. 2965 

5510 

2. 91 X 10' 5 

2535 

0. 867 

11320 

3. 08 X 10" 5 



Ave rage : 

3. 10 X 10~ 5 
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TABLE IV 


TYPICAL ELECTRON TEMPERATURE MEASUREMENTS 


AS A FUNCTION OF RF POWER FOR NITROGEN 

RF Power: 0 watts 

4. 2 w 

8. 4 w 

12. 6 w 

16. 8 w 

Electron Temperature, °K 

3500 

4320 




3490 

4790 




3390 

45™ 





4590 

5950 




4640 

5710 





6110 

7510 


3540 



7340 


3410 




8940 

3550 




8680 

Averages: 3480 

4580 

5920 

7430 

8810 
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TABLE V 


TYPICAL ENERGY-LOSS FACTORS FOR NITROGEN 
AS A FUNCTION OF GAS TEMPERATURE 


o T . 

Gas Ten perature, K 
V T ^ T e,E = 0 

rf Power, 
watts 

Electron Temp, , 
T g , °K (with rf ) 

Energy-loss 
Factor, 6 gff 

1260 

2. 32 

3865 

2. 54 X 10‘ 4 

1700 

2. 4 

4130 

2. 9 X 10‘ 4 

1715 

2. 24 

3915 

2. 92 X 10" 4 

1750 (Air leak) 

4. 2 

3260 

7.8 X 10‘ 4 

2860 

2. 4 

3950 

6. 3 X 10" 4 

2940 

2. 4 

3910 

7. 1 X 10" 4 

3190 

2. 4 

4315 

6. 1 X 10” 4 

3190 

0. 8 

3500 

7. 6 X 10“ 4 

3190 

2. 4 

4250 

6. 5 X 10' 4 

3480 (Air leak) 

4. 2 

4580 

10. 9 X 10' 4 

3610 

2. 4 

4600 

6. 95 X 10' 4 

3760 

2. 4 

4790 

6.7 X 10‘ 4 

3830 

2. 24 

4650 

7. 9 X 10“ 4 

3970 

0. 8 

4320 

6.5 X 10' 4 

4280 

2. 4 

5340 

6.7 X 10' 4 

5026 

2. 4 

6065 

6. 6 X 10" 4 

5950 

2. 4 

6490 

13. 0 X 10* 4 

6100 

2. 4 

7540 

4. 9 X 10 -4 
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TABLE VI 


TYPICAL ENERGY-LOSS FACTORS FOR OXYGEN- NITROGEN MIXTURES AS A FUNCTION OF 

COMPOSITION AND CAS AND ELECTRON TEMPERATURES (R~ = Q .. /Q _ ) 

Q e ^2 e< - > 2 


Derived energy loss factor 


[o 2 ] 

[o 2 ] +[N 2 ] 

Gas temp. °K 

T = T 
g e f E = 0 

rf power 
W 

Mixture 

Electron temp. energy-loss 

T ,°K (with rf) factor, & e ^X 10 

for pure , 

using 6 n = 

4 ^ 

R q = 2.5 

6 ° 2 x !° 
3 X 10 

R q = 2. 

0 

4060 

2. 4 

4920 

8. 0 



0 

3240 

2. 4 

3950 

9. 7 



0 

2960 

4. 8 

4700 

7. 92 



0 

3530 

1. 6 

4280 

6. 12 



0. 0356 

2220 

2. 4 

3210 

6. 96 

274 

221 

0. 0416 

2030 

2. 4 

3160 

6. 10 

184 

148 

0. 0453 

1750 

4. 8 

3415 

8. 27 

286 

230 

0. 0453 

1815 

8. 0 

3870 

11. 2 

442 

356 

0. 0834 

2000 

2.4 

2550 

12. 5 

274 

222 

0. 0834 

1730 

4. 8 

2880 

12. 0 

256 

208 

0. 0834 

1850 

8. 0 

3360 

15. 2 

346 

284 

0. 1007 

1630 

4. 8 

2550 

15. 0 

286 

231 

0. 1251 

1680 

2. * 

1920 

28. 7 

476 

389 

0. 1431 

1850 

8. 0 

2980 

20. 3 

279 

228 

0. 1666 

1690 

2. 4 

1970 

24. 6 

294 

241 

0. 1666 

1865 

8. 0 

2520 

35. 0 

435 

355 



TABLE VI 
(Concluded) 


0. 1850 

1240 

2. 4 

1580 

20. 3 

211 

173 

0. 2217 

1680 

8. 0 

2100 

54. 7 

506 

417 

0. 2217 

1730 

16. 0 

2440 

64. 7 

606 

498 

0. 2481 

1 660 

2. 4 

1820 

43. 0 

346 

285 

0. 2481 

1430 

2. 4 

1700 

25. 5 

196 

162 

0. 3006 

1440 

2. 4 

1810 

18. 6 

109 

91 

0. 3211 

1540 

16. 0 

2520 

46. 9 

279 

233 

0. 3211 

1630 

28. 0 

2670 

77. 3 

470 

391 

0. 3333 

1590 

2.4 

1750 

43. 0 

243 

203 

0. 4169 

1495 

2. 4 

1700 

33. 6 

145 

119 

0. 500 

1275 

16 

2050 

59. 3 

200 

173 

0. 500 

1760 

56 

4550 

57. 6 

194 

167 

0. 500 

1595 

2. 4 

1865 

25. 5 

82 

71 





Average : 

297 

244 



TABLE VII 


DEPENDENCE OF ENERGY LOSS FACTOR FOR PURE OXYGEN, 
DERIVED FROM MIXTURE DATA BY MIXTURE RULE 1 AND 
ASSUMING THAT ONLY N 2 AND 0 2 ARE PRESENT IN THE 
MIXTURE, ON THE OXYGEN CONTENT OF THE MIXTURE (R Q = 2. 5) 


[o 2 ] 

4 

Derived 6 n X 10 
U 2 

[o 2 l 

Mn 2 ] 

~ 0. 04 

(four runs) 

296 

0. 083 

(three runs) 

292 

~0. 16 

(four runs) 

305 

~0. 23 

(four runs) 

413 

~0. 32 

(four runs) 

275 

0. 50 

(three runs) 

159 
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2. 45 GHz Illumination System 
Item Description 

1 rf Power Source Raytheon J 

PGM-10-X-2 (85 Watts f 

Max) 

2 Directional Coupler - NARDA 

3003-30 

3 Stub Tuner - NARDA 903 N 

4 Load -NARDA 369 NM 

5 20 DB PAD-NARDA 757-20 

6 Thermistor Mount HP 447 B 

7 Power Meter H. P. 430 C 

8 Vacuum Tank 

9 Waveguide Test Section 

10 Variable Attenuator 



Figure 1. - Schematic of microwave plasma illumination system. 






8 eff For N 2 - 0 2 Mixtures ot 
Tg Ty® Xe t E - 0^ I700°K ond 
T e between !500 # Kond 4500°K- 


Pure N 2 C99.5°/< 
at Tg « 3450° K 


mm 


Pure N 2 at Tg » T v » T e ,E-o < *800° K 


0.04 


Figure 3, Energy-loss factor for mixtures at 

T « T « T « 1700 °K as a function 

g vib e, E = 0 

of oxygen mole fraction. Purity of is 

99. 5% as claimed by manufacturer. 






Figure 1, Appendix C. Langmuir probe characteristics for nitrogen 

plasmajet with 2. 4 watts and 0 rf power 
(traced twice). 
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Figure 3, Appendix C. Radial survey of electron and ion current to 

Langmuir probe obtained by continuous traverses 
of probe through plasmajet axis. 





Electron current 



0.4 Axis 


0.01 


Ion current 


R = 100 Ohm 


Figure 4, Appendix C. Axial survey of electron and ion current to Langmuir probe 

obtained by continuous traverses of probe along plasmajet 
axis. Peaks in electron current correspond to passage of 
probe through shock diamonds. 
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Figure 6, Appendix C. 


obtained 


r probe characteristics along axis of plasmajet 
at l" intervals from nozzle exit to 20" downstream. 










Figure 8, Appendix C. Effect of RF on probe electron current on jet axis (axial 

traverse). Note the effect is negligible outside waveguide 
(two traces come together). 
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Figure 9, Appendix C. Rotational band of N 
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rotational temperature. 
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